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Following the general framework of the gauge invariant perturbation theory developed in 
the papers [K. Nakamura, Prog. Theor. Phys. 110 (2003), 723; ihid. 113 (2005), 481], we for- 
mulate second-order gauge invariant cosmological perturbation theory in a four- dimensional 
homogeneous isotropic universe. We consider perturbations both in the universe dominated 
by a single perfect fluid and in that dominated by a single scalar field. We derive all the 
components of the Einstein equations in the case that the first-order vector and tensor modes 
are negligible. All equations are derived in terms of gauge invariant variables without any 
gauge fixing. These equations imply that second-order vector and tensor modes may be 
generated due to the mode-mode coupling of the linear-order scalar perturbations. We also 
briefiy discuss the main progress of this work through comparison with previous works. 

§1. Introduction 

The general relativistic cosmological linear perturbation theory has been devel- 
oped to a high degree of sophistication during the last 25 years. One of the 
motivations of this development is to clarify the relation between the scenarios of 
the early universe and cosmological data, such as the cosmic microwave background 
(CMB) anisotropics. Recently, the first-order approximation of our universe from 
a homogeneous isotropic one was revealed through the observation of the CMB by 
the Wilkinson Microwave Anisotropy Probe (WMAP).^^ This observation suggests 
that the fluctuations are adiabatic and Gaussian at least to a first order approxi- 
mation. One of the next important theoretical studies is to clarify the accuracy of 
these results by studying the non-Gaussian behavior, non-adiabatic behaviors, and 
so on. These will be goals of future satellite missions. To estimate the accuracy 
of the first-order approximation, theoretically, it is necessary to investigate second- 
order cosmological perturbations. From the observational point of view, also, with 
the increase of precision of the CMB data, the study of relativistic cosmological per- 
turbations beyond linear order is becoming a topical subject, especially in regard 
to studying the generation of the primordial non-Gaussian behavior in inflationary 
scenarios^^ and the non-Gaussian component in the CMB anisotropy. 

In the literature, the second-order general relativistic perturbation theory has 
been investigated by many researchers. For the pioneering work, Tomita^^ investi- 
gated general relativistic second-order perturbations in the Einstein-de Sitter model 
(vanishing A model), and his treatment is in the synchronous gauge. His second- 
order perturbation theory was later extended to the general relativistic Zel'dovich ap- 
proximation. ^-^ Recently, nonlinear gauge transformations and the concept of gauge 
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invariance have been studied by Bruni et al.^-* As the by-product of their research, 
Sonego and Bruni^^^ obtained a representation of the higher-order Taylor expansion 
of tensors on a manifold in a quite generic form, and the second-order gauge trans- 
formation from the synchronous gauge to the Poisson gauge has been performed 
by Matarrese et al.^^^ More recently, Noh and Hwang^^^ derived second-order per- 
turbation equations in the Priedmann universe. Further, Tomita^^^ also extended 
his original works to second-order perturbations of nonzero- /I cosmological mod- 
els and studied the CMB anisotropy and its non-Gaussian nature. However, their 
treatments of the general relativistic second-order perturbation is very complicated. 
Hence, to avoid this complicated formulation, there have also been several attempts 
to investigate the nonlinear effects of general relativistic perturbations. 

In this paper, we present a very clear formulation of the general relativistic 
second-order cosmological perturbations in a homogeneous isotropic universe. This 
paper is the complete version of the previous short paper^^^ by the present author. 
Our formulation in this paper is one of the applications of the gauge invariant formu- 
lation of the second-order perturbation theory on a generic background spacetime 
developed in two papers by the present author. These papers are referred to in 
this paper as KN2003^^) and KN2005.^'^) This formulation is a by-product of inves- 
tigations of the oscillatory behavior of a self-gravitating Nambu-Goto membrane^^^ 
and was first applied to a comparison between the oscillatory behavior of a gravitat- 
ing Nambu-Goto string and that of a test string.^^^ This was a trivial application of 
the general formulation developed in KN2003 and KN2005, while the second-order 
cosmological perturbation carried out in this paper is the first non-trivial application 
of the formulation presented in KN2003 and KN2005. 

The formulation developed in KN2003 and KN2005 is an extension of the works 
of Bruni et al.^^ The gauge transformation rules of the perturbations formulated by 
Bruni et al.^-* are extended to those in multi-parameter perturbation theory. ^^^'^'^^ 
Based on these gauge transformation rules, in KN2003, we proposed a procedure 
to find gauge invariant variables on a generic background spacetime to third-order 
perturbations, assuming that we already know the procedure to find gauge invariant 
variables for the linear-order metric perturbations. We also showed in KN2005 that 
the proposal of the gauge invariant variables in KN2003 provides a self-consistent 
second order perturbation theory in the generic background spacetime. It is straight- 
forward to apply this general formulation to cosmological perturbations. In the 
cosmological perturbation case, there are some proposals of gauge invariant formu- 
lations of the second-order perturbation. For example, Mukhanov et al.^^^ proposed 
a gauge invariant second-order perturbation to evaluate the back reaction effect of 
the inhomogeneities in the universe on the effective expansion law of the universe. 
However, we should distinguish our formulation presented in this paper from the 
proposal of Mukhanov et al.^^^ These are quite different approaches. 

To develop the gauge invariant perturbation theory, we start by explaining the 
concept of the "gauge" in general relativistic perturbation theory to avoid any mis- 
understanding of our formulation. General relativity is based on the concept of 
general covariance. Intuitively, the principle of general covariance states that there 
is no preferred coordinate system in nature, though the notion of general covariance 
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is mathematically included in the definition of a spacetime manifold in a trivial way. 
This is based on the philosophy that coordinate systems are originally chosen by us, 
and natural phenomena have nothing to do with our coordinate systems. If we apply 
a peculiar coordinate system to investigate natural phenomena, we will see peculiar 
behavior of that natural phenomena due to this peculiar coordinate system. This is 
an intuitive explanation of general covariance. Due to this general covariance, the 
gauge degree of freedom, which is the unphysical degree of freedom of perturbations, 
arises in general relativistic perturbations. To obtain physically meaningful results, 
we have to fix this gauge degree of freedom or to extract the gauge invariant part of 
perturbations. 

As reviewed in detail in ^ of this paper, the developments in KN2003 and 
KN2005 are based on the understanding of the "gauge" in the perturbation theory 
which was first proposed by Stewart et al.^^^ and developed by Bruni et al.^^'^'''' 
Based on this formulation, we define the complete set of gauge invariant variables 
of the second-order cosmological perturbations in the Friedmann-Robertson- Walker 
universe. We consider two cases of the Friedmann-Robertson- Walker universe, one 
in which the universe is filled with a single perfect fiuid and one in which the universe 
is filled with a single scalar field. We also derive the second-order Einstein equations 
of cosmological perturbations in terms of these gauge invariant variables without any 
gauge fixing in these two cases. 

The organization of this paper is as follows. In fj2l we review the general frame- 
work of the second-order gauge invariant perturbation theory developed in KN2003 
and KN2005. This review also includes additional explanation not given in those 
papers. In ij3l we summarize the Einstein equations in the case of a background 
homogeneous isotropic universe, which are used in the derivation of the first- and 
second-order Einstein equations. In ^ we define the first- and second-order gauge 
invariant variables for the cosmological perturbations. The first-order perturbation 
of the Einstein equations is reviewed in ^ This perturbation is used in the deriva- 
tion of the second-order Einstein equations. Then, the derivation of the second-order 
Einstein equation is given in ^ The final section, fJTl is devoted to a summary and 
discussions concerning the relation between this and previous works. 

We employ the notation of KN2003 and KN2005 and use abstract index nota- 
tion. We also employ natural units in which Newton's gravitational constant is 
denoted by G and the velocity of light satisfies c = 1. 

§2. General framework of the gauge invariant perturbation theory 

In this section, we briefiy review the general framework of the gauge invariant 
perturbation theory developed in KN2003 and KN2005 by the present author. To 
explain the gauge degree of freedom in perturbation theories, we have to recall what 
we are doing when we consider perturbations. Further, we comment on the Taylor 
expansion of tensors on a manifold, at first, in ^2.H because any perturbation theory 
is based on the Taylor expansion. Next, in ^2.21 we review the basic understanding of 
the gauge degree of freedom in perturbation theory based on the work of Stewart et 
al.^^^ and Bruni et al.^^ When we consider perturbations in the theory with general 
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covariance, we have to exclude these gauge degrees of freedom in the perturbations. 
To accomphsh this, gauge invariant variables of perturbations are useful, and these 
are regarded as physically meaningful quantities. In ^2.3^ we review the procedure 
for finding gauge invariant variables of perturbations, which was developed by the 
present author in KN2003. Then, in ^2.41 we briefly review the general issue of 
the gauge invariant formulation for the second-order perturbation of the Einstein 
equation developed in KN2005. We emphasize that the ingredients of this section 
do not depend on the background spacetime, and they are applicable not only to 
cosmological perturbations but also to any other general relativistic perturbations. 

2.1. Taylor expansion of tensors on a manifold 

Here, we comment on the general form of the Taylor expansion of tensors on a 
manifold A4. We first consider the Taylor expansion of a scalar function / : i— > M, 
which can be extended to any tensor field on a manifold. 

The Taylor expansion of a function / is an approximated form of f{q) at g G 
in terms of the variables at p G Ai, where q is in the neighborhood of p. To 
consider the Taylor expansion of a function /, we introduce a one-parameter family 
of diffeomorphisms (p\ : A4 i— > 7W, where ^x{p) = q and ^x=o{p) = p. One example 
of a diffeomorphisms (Px is an exponential map. However, we consider a more general 
class of diffeomorphisms. 

In terms of the diffeomorphism ^x^ the Taylor expansion of the function f[q) is 
given by 
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A2+0(A3). (2-1) 



Since this expression hold for an arbitrary smooth function /, we may regard the 
Taylor expansion to be the expansion of the pull-back (p*^ of the diffeomorphism ^x , 
rather than the expansion of the function /. 

Further, as shown by Sonego and Bruni,^'^^ there exist vector fields and ^2 
such that the expansion ()2Tp is given by 

f{q) = {^lf){p) = f{p) + {£U\ A + i {£^, + £l) fl A2 + 0(A3), (2-2) 

without loss of generality. In the representation (l2-2p of the Taylor expansion, 
and ^2 the generators of the one-parameter family of diffeomorphisms ^x and 
these represent the direction along which the Taylor expansion is carried out. The 
generator is the first-order approximation of the fiow of the diffeomorphism <Px^ 
and the generator ^2 is the second-order correction to this flow. 

When the generator ^2 is proportional to the generator the representation 
(|2-2p of the Taylor expansion reduces to that of the pull-back of an exponential map. 
Therefore, we may regard that the Taylor expansion (j2-2p is the generalization of 
an exponential map (one-parameter group of diffeomorphisms). However, as shown 
by Sonego and Bruni,^'^^ the Taylor expansion of a C" one-parameter family of 
diffeomorphisms can always be represented in the form of ()2-2p . In general, the 
generator ^2 may not be proportional to the generator ^J. Hence, we regard the 
generators and to be independent. 
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Further, ^\ can be extended to diffeomorphisms acting on tensor fields of all 
types. Thus, the Taylor expansion of a tensor field Q of any type on a manifold M. 
is given by 

Qiq) = {mm = QiP) + A + ^ + £l) Q^X^ + 0(A3), (2-3) 

and we conclude that the representation of the Taylor expansion (j2-3p is quite general. 

2.2. Gauge degree of freedom in perturbation theory 

Now, we explain the concept of gauge in general relativistic perturbation theory. 
To explain this, we first point out that, in any perturbation theory, we always treat 
two spacetime manifolds. One is the physical spacetime A4, which we attempt to 
describe in forms of perturbations, and the other is the background spacetime Mq, 
which is a fictitious manifold which we prepare for perturbative analyses. We empha- 
size that these two spacetime manifolds M and Mq are distinct. Let us denote the 
physical spacetime by {M.,gab) and the background spacetime by (Moidab), where 
(jab is the metric on the physical spacetime manifold, Ai, and gab is the metric on 
the background spacetime manifold, A4q. Further, we formally denote the spacetime 
metric and the other physical tensor fields on the physical spacetime by Q and its 
background value on the background spacetime by Qq. 

Second, in any perturbation theories, we always write equations for the pertur- 
bation of the physical variable Q in the form 

QVy) = Qo{p) + SQ{p). (2-4) 

Usually, this equation is simply regarded as a relation between the physical variable 
Q and its background value Qo, or as the definition of the deviation 6Q of the 
physical variable Q from its background value Qq. However, Eq. (j2-4p has deeper 
implications. Keeping in our mind that we always treat two different spacetimes, 
{^A,gab) and {M.o,gab), in perturbation theory, Eq. (j2-4p is a rather strange equation 
in the following sense: The variable on the left-hand side of Eq. (j2-4p is a variable 
on the physical spacetime {M,gab), while the variables on the right-hand side of 
Eq. (|2-4p are varibles on the background spacetime, {Mo, dab)- Hence, Eq. (j2-4p 
gives a relation between variables on two different manifolds. 

Further, through Eq. (j2-4p . we have implicitly identified points in these two 
different manifolds. More specifically, QC'p") on the left-hand side of Eq. (j2-4p is a 
field on the physical spacetime, A4, and "p" G A4. Similarly, we should regard the 
background value Qo{p) of Q{^^p^^) and its deviation 5Q{p) of from Qq{p), 

which are on the right-hand side of Eq. (j2-4p . as fields on the background spacetime, 
A^O) and p G Mq. Because Eq. ()2-4p is regarded as an equation for field variables, 
it implicitly states that the points "p" € M. and p € TWq are same. This represents 
the implicit assumption of the existence of a map A^o M. : p & M.Q ^ "p" G A^, 
which is usually called a gauge choice in perturbation theory. ^^-^ 

It is important to note that the correspondence between points on Mq and A4, 
which is established by such a relation as Eq. (|2-4p . is not unique to the theory 
in which general covariance is imposed. Rather, Eq. ()2-4p involves the degree of 
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freedom corresponding to the choice of the map X : A4q A4. This is called the 
gauge degree of freedom. Such a degree of freedom always exists in perturbations of 
a theory in which we impose general covariance. "General covariance" intuitively 
means that there is no preferred coordinate system in the theory. If general covari- 
ance is not imposed on the theory, there is a preferred coordinate system, and we 
naturally introduce this coordinate system onto both A^o and M. Then, we can 
choose the identification map X using this coordinate system. However, there is 
no such coordinate system in general relativity, due to its general covariance, and 
we have no guiding principle to choose the identification map X. Indeed, we could 
identify "p" G A4 with q £ M.q {q ^ p) instead of p € Mq. In the above understand- 
ing of the concept of "gauge" in general relativistic perturbation theory, a gauge 
transformation is simply a change of the map X. 

These are the basic ideas necessary to understand gauge degree of freedom in the 
general relativistic perturbation theory proposed by Stewart and Walker. This 
understanding has been developed by Bruni et al.,^-* and by the present author. 
We briefly review this development. 

To formulate the above understanding in more detail, we introduce an infinites- 
imal parameter A for the perturbation. Further, we consider the 4 + 1-dimensional 
manifold = x R, where 4 = dimAi and A G R. The background spacetime 
A^o = -^Ia=o physical spacetime M. = M\ = N\^^^^ are also submanifolds 

embedded in the extended manifold N . Each point on J\f is identified by a pair, 
(p. A), where p € AIa, and each point in the background spacetime Aio in is 
identified by A = 0. 

Through this construction, the manifold TV is foliated by four-dimensional sub- 
manifolds AIa of each A, and these are diffeomorphic to the physical spacetime Ai 
and the background spacetime A^o- The manifold TV has a natural differentiable 
structure consisting of the direct product of TW and R. Further, the perturbed 
spacetimes M.\ for each A must have the same differential structure with this con- 
struction. In other words, we require that perturbations be continuous in the sense 
that (TW, and (A^Oi 5afe) are connected by a continuous curve within the extended 
manifold TV. Hence, the changes of the differential structure resulting from the per- 
turbation, for example the formation of singularities and singular perturbations in 
the sense of fluid mechanics, are excluded from consideration in this paper. 

Let us consider the set of field equations 

£[Qx] = (2-5) 

on the physical spacetime M.\ for the physical variables Q\ on A4\. The field 
equation (|2-5|) formally represents the Einstein equation for the metric on ^A\ and 
the equations for matter fields on If a tensor field Q\ is given on each 

Qx is automatically extended to a tensor field on TV by Q{p,X) := Qx{p), where 
p € A4x. In this extension, the field equation (|2-5|) is regarded as an equation on 
the extended manifold TV. Thus, we have extended an arbitrary tensor field and the 
field equations (j2-5p on each A4x to those on the extended manifold TV. 

Tensor fields on TV obtained through the above construction are necessarily "tan- 
gent" to each A4x, i.e., their normal component to each Mx identically vanishes. To 
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consider the basis of the tangent space of we introduce the normal form and its 
dual, which are normal to each M.\ in M. These are denoted by {dX)a and ((9/9A)", 
respectively, and they satisfy 

(^^)<^ fIrV = 1- (2-6) 



The form {d\)a and its dual, are normal to any tensor field extended from 

the tangent space on each A4\ through the above construction. The set consisting 
of {d\)a-, {d/dXY and the basis of the tangent space on each M.\ is regarded as the 
basis of the tangent space of M. 

To define the perturbation of an arbitrary tensor field Q, we compare Q on the 
physical spacetime Aix with Qo on the background spacetime, and it is necessary 
to identify the points of M.\ with those of Mq. This point identification map is the 
so-called gauge choice in the context of perturbation theories, as mentioned above. 
The gauge choice is made by assigning a diffeomorphism X\ : J\f ^ M such that X\ 
: Aio M-x- Following the paper of Bruni et al.,^) we introduce a gauge choice Xx 
as one of the one-parameter groups of diffeomorphisms, i.e., an exponential map, for 
simplicity. We denote the generator of this exponential map by "^r;". This generator 
■^77" is decomposed by the basis on the tangent space of Af which are constructed 
above. Though the generator '^rj"' should satisfy some appropriate conditions, the 
arbitrariness of the gauge choice Xx is represented by the tangential component of 
the generator '^rj"' to the tangent space of A4x- 

The pull-back X^Q, which is induced by the exponential map Xx, maps a tensor 
field Q on the physical manifold A4x to a, tensor field X^Q on the background 
spacetime. In terms of this generator '^rj"', the pull-back X^Q is represented by the 
Taylor expansion 

Q{r) = Q{Xx{v)) = X*xQ{p) = Q(P) + A £^r,Ql + £1^q\^ + 0(A3), (2-7) 
where r = Xx{p) Mx- Because p S Mo^ we may regard the equation 

X*xQ{p) = Qoip) + A £^,Q\j^^ ip) + ^a2 £%Q\^^ (p) + O(A^) (2-8) 

as an equation on the background spacetime A^O) where Qq = Q\j^^^ is the back- 
ground value of the physical variable of Q. Once the definition of the pull-back of 
the gauge choice Xx is given, the perturbation A^Qx of a tensor field Q under the 
gauge choice Xx is simply defined as 

A^Qx := ^A*QUo - Qo- (2-9) 

We note that all variables in this definition are defined on A4q. Expanding the first 
term on the right-hand side of (j2-9p as 



^xQxImo =Qo + A^i^Q + ^A^^Jq + 0(A3), (2-10) 
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we define the first- and the second-order perturbations of a physical variable Q\ 
under the gauge choice Xx by 



(2-11) 



Now, we consider two different gauge choices based on the above understanding 
of the gauge choice in the perturbation theory. Suppose that X\ and y\ are two 
exponential maps with the generators "^rf and ^rf on A/", respectively. In other 
words, X\ and y\ are two gauge choices. Then, the integral curves of each '^•rf and 
^jf in are the orbits of the actions of the gauge choices X\ and ^^a, respectively. 
Since we choose the generators '^'rf and ^'rf so that these are transverse to each 
}A\ everywhere on AA, the integral curves of these vector fields intersect with each 
Therefore, points lying on the same integral curve of either of the two are 
to be regarded as the same point within the respective gauges. When these curves 
are not identical, i.e., the tangential components to each Mx of "^r/" and ■^r/'* are 
different, these point identification maps Xx and yx are regarded as two different 
gauge choices. 

We next introduce the concept of gauge invariance. Following the paper by 
Bruni et al.,^''^ we consider the concept of gauge invariance up to order n. Suppose 
that Xx and yx are two different gauge choices which are generated by the vector 
fields ^r]"" and ^7]°", respectively. These gauge choices also pull back a generic tensor 
field Q on A/" to two other tensor fields, X^Q and J^^Q, for any given value of A. In 
particular, on A^O; we now have three tensor fields associated with a tensor field Q; 
one is the background value Qq of Q, and the other two are the pulled-back variables 
of Q from M.x to A4o by the two different gauge choices, 

■^Qa := -^a*QIa4o = Q^ + >'^xQ + \>^^^xQ + 0{X') =Qo + A^Qx. (2-12) 

^Qa := yiQ\M, = + A^i^Q + \\^^§Q + 0{>?) = Qo + ^^Qa- (2-13) 

Here, we have used Eqs. (|2-9p and (j2-10p . Because Xx and yx are gauge choices that 
map the background spacetime M.q into the physical spacetime A4a, ^Qx and ^Qx 
are the representations on A^o of the perturbed tensor field Q in the two different 
gauges. The quantities and '^yQ in Eqs. (|2-12|) and (|2-13|) are the perturbations 
of 0{k) in the gauges Xx and ^^a, respectively. We say that Q is gauge invariant up 
to order n iff for any two gauges Xx and 3^a the following holds: 

(J)q = ^'^Q yk, with k<n. (2-14) 

From this definition, we can prove that the nth-order perturbation of a tensor field 

(k) 

Q is gauge invariant up to order n iff in a given gauge Xx we have £^\^Q = for 
any vector field ^'^ defined on Mq and for any k < n. As a consequence, the nth- 
order perturbation of a tensor field Q is gauge invariant up to order n iff Qq and all 
its perturbations of lower than nth order are, in any gauge, either vanishing or con- 
stant scalars, or a combination of Kronecker deltas with constant coefficients.^^' ^''-''^^^ 
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Further, even if its lower order perturbations are not trivial, we can decompose any 
perturbation of Q into the gauge invariant and gauge variant parts, as shown in 
KN2003. This will be also explained in the next subsection. 

Now, we consider the gauge transformation rules between different gauge choices. 
In general, the representation "^Qx on M.q of the perturbed variable Q on M.x de- 
pends on the gauge choice X\. If we employ a different gauge choice, the representa- 
tion of Q\ on M.0 may change. Suppose that X\ and y\ are different gauge choices, 
which are the point identification maps from Mq to M.x, and the generators of these 
gauge choices are given by "^r/" and -^ry", respectively. Then, the change of the gauge 
choice from X\ to y\ is represented by the diffeomorphism 

c!>^:={Xx)-^oy,^. (2-15) 

This diffeomorphism ^\ is the map ^\ : A4q — > A4q for each value of A G M. The 
diffeomorphism does change the point identification, as expected from the under- 
standing of the gauge choice discussed above. Therefore, the diffeomorphism (Px is 
regarded as the gauge transformation (Px '■ Xx yx- 

The gauge transformation $x induces a pull-back from the representation "^Qx 
of the perturbed tensor field Q in the gauge choice Xx to the representation ^Qx in 
the gauge choice yx- Actually, the tensor fields "^Qx and ^Qxj which are defined on 
TWo; are connected by the linear map (p\ as 



Mo 



= {x^'yxYixiQ) =<^l''Qx- (2-16) 

According to generic arguments concerning the Taylor expansion of the pull-back of 
a tensor field on the same manifold, given in ^2.11 it should be possible to express 
the gauge transformation 'P\'^Qx in the form 

^l^Q = ^Q + A%^Q + y + 4J + 0(A3), (2-17) 

where the vector fields £,f and ^2 generators of the gauge transformation <Px- 

Comparing the representation (j2-17p of the Taylor expansion in terms of the 
generators and .^f of the pull-back (^^"^Q and that in terms of the generators 
■^r/" and ^77" of the pull-back y^ o {X^^)* ■^Q (= (P*x'^Q), we readily obtain explicit 
expressions for the generators and of the gauge transformation ^ = X^"^ o y^^ 
in terms of the generators '^r/" and -^r/" of the gauge choices as follows: 

e? = V-V, = [y^,•^r/]^ (2-18) 

Further, because the gauge transformation <Px is a map within the background space- 
time Alo) the generator should consist of vector fields on Mq. This can be satisfied 
by imposing some appropriate conditions on the generators •^r/'^ and '^7]°'. 

We can now derive the relation between the perturbations in the two different 
gauges. Up to second order, these relations are derived by substituting (|2-12p and 
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(1243]) into (12^71) : 



2-^5(1) ^xQ + { ■^5(2) + "^l(i) } 




'0- 



(2-19) 
(2-20) 



These results are, of course, consistent with the concept of gauge invariance up to 
order n, as introduced above. Further, inspecting these gauge transformation rules, 
we can define the gauge invariant variables at each order, as shown below. 

Here, we should comment on the gauge choice in the above explanation. We 
have introduced an exponential map X\ (or yx) as the gauge choice, for simplicity. 
However, this simplified introduction of the gauge choice X\ as an exponential map 
is not essential to the gauge transformation rules (|2-19p and (|2-20|) . Indeed, we 
can generalize the diffeomorphism X\ from an exponential map. If we generalize 
the diffeomorphism X\, the representation (j2-8p of the pulled-back variable X^Q{p), 
the representations of the perturbations (|2-lip . and the relations (|2-18p between 
generators of Xx, and yx will be changed. However, the gauge transformation 
rules (|2-19p and (j2-20p are direct consequences of the Taylor expansion (j2T7p of <Px ■ 
As commented in ^2.11 the representation of the Taylor expansion (j2-17p of ^x is 
quite general. ^'^^ Therefore, the gauge transformation rules (j2T9p and (j2-20p do not 
change, even if we generalize the choice of Xx- Further, the relations (|2-18p between 
generators also imply that, even if we consider a simple exponential map as the gauge 
choice, both of the generators and ^2 ^■re naturally induced by the generators of 
the original gauge choices. Hence, we conclude that the gauge transformation rules 
(|2-19p and (j2-20p are quite general and irreducible. In this paper, we develop a 
second-order gauge invariant cosmological perturbation theory based on the above 
understanding of the gauge degree of freedom only through the gauge transformation 
rules (j2-19p and (|2-20p . Hence, the development of the cosmological perturbation 
theory presented below is not changed if we generalize the gauge choice Xx from a 
simple exponential map. 

2.3. Gauge invariant variables 

Inspecting the gauge transformation rules (j2-19p and (j2-20p . we can define the 
gauge invariant variables for a metric perturbation and for arbitrary matter fields. 
Employing the idea of gauge invariance up to order n for nth-order perturbations,^^ 
we proposed a procedure to construct gauge invariant variables of higher-order per- 
turbations. This proposal is as follows. First, we decompose a linear-order metric 
perturbation into its gauge invariant and variant parts. The procedure for decom- 
posing linear-order metric perturbations is extended to second-order metric pertur- 
bations, and we can decompose the second-order metric perturbations through a 
procedure similar to that for the linear-order metric perturbation. Then, we de- 
fine the gauge invariant variables for the first- and second-order perturbations of an 
arbitrary field other than the metric by using the gauge variant parts of the first- 
and second-order metric perturbations. Though the procedure for finding gauge 
invariant variables for linear-order metric perturbations is highly non-trivial, once 
we know this procedure, we can easily find the gauge invariant part of a higher- 
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order perturbation through a simple extension of the procedure for the hnear-order 
perturbations. 

To consider a metric perturbation, we expand the metric on the physical space- 
time M., which is pulled back to the background spacetime A^o using a gauge choice 
in the form given in (j2-10p . 

^X9ab = 9ab + Xxhab + y Af, + 0%X), (2-21) 

where gab is the metric on the background spacetime A4q. Of course, the expansion 
(l2-2ip of the metric depends entirely on the gauge choice X\. Nevertheless, hence- 
forth, we do not explicitly express the index of the gauge choice X\ in an expression 
if there is no possibility of confusion. 

Our starting point to construct gauge invariant variables is the assumption that 
we already know the procedure for finding gauge invariant variables for the linear 
metric perturbations. Then, a linear metric perturbation hab is decomposed as 

hab =■ 'Hab + £xgab, (2-22) 

where Tiab and X" are the gauge invariant and variant parts of the linear-order metric 
perturbations,^^-* i.e., under the gauge transformation p-19p . these are transformed 
as 

yHab - xHab = 0, yX°- - xX"" = i'^^y (2-23) 

As emphasized in KN2003 and KN2005, the above assumption is quite strong 
and it is not simple to carry out the systematic decomposition (j2-22p on an arbitrary 
background spacetime, since this procedure depends completely on the background 
spacetime (A^o,5ab)- However, we show that this procedure exists in the case of 
cosmological perturbations of a homogeneous and isotropic universe in ^A.W 

Once we accept this assumption for linear-order metric perturbations, we can 
always find gauge invariant variables for higher-order perturbations. ^^-^ As shown in 
KN2003, at second order, the metric perturbations are decomposed as 

hb ='■ ^ab + '^£xhab + {£y - £x) gab, (2-24) 

where Cab and Y"" are the gauge invariant and variant parts of the second order 
metric perturbations, i.e., 

yJCab - xCab = 0, 3;r"-A'l^" = e(2) + [e(i),^r. (2-25) 

The details of the derivation of this gauge invariant part of the second-order metric 
perturbation are explained in the context of cosmological perturbations in ^4.2[ 

Furthermore, as shown in KN2003, using the first- and second-order gauge vari- 
ant parts, X"" and y", of the metric perturbations, the gauge invariant variables for 
an arbitrary field Q other than the metric are given by 

Wq := - £^Qo, (2-26) 
:= (2)Q _ 2£x^^'^Q - {£y - £x} Qo- (2-27) 
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It is straightforward to confirm that the variables ^^^Q defined by (j2-26p and (j2-27p are 
gauge invariant under the gauge transformation rules (j2-19p and (j2-20p . respectively. 

Equations (j2-26p and (j2-27p have several more important implications. To see 
this, we represent these equations as 

(1) q = + £xQo, (2-28) 

(2) q = (2)q + 2£x^^^Q + {£y- £x}Qo- (2-29) 

These equations imply that any perturbation of first and second order can always be 
decomposed into gauge invariant and gauge variant parts as Eqs. (|2-28p and (|2-29p . 
respectively. In ^4.31 we see that these formulae for the decomposition into gauge 
invariant and variant parts of each order perturbation are very important. 

2.4. Perturbations of the Einstein tensor and the Einstein equations 

Now, we review the formulae for the perturbative Einstein tensor at each order 
that are presented in KN2005. The relation between the curvatures associated with 
the metrics on the physical spacetime Aix and the background spacetime A4o is 
given by the relation between the pulled-back operator X^Va {X\^) of the covari- 
ant derivative associated with the metric Qab on Mx and the covariant derivative 
Vq associated with the metric gab on Mq. The pulled-back covariant derivative 
A'j^Va {X\^) depends on the gauge choice Xx. The property of the derivative oper- 
ator X^Va {Xx^Y as the covariant derivative on the physical spacetime Mx is given 
by 

X^Va [{X-Yx^9ab) = 0, (2-30) 

where X^gab is the pull-back of the metric on the physical spacetime Aix, which is ex- 
panded as Eq. (j2-2ip . In spite of the gauge dependence of the operator X^Va (X^^)* , 
we simply denote this operator by Va, because our calculations are carried out only 
on the background spacetime A^o in the same gauge choice Xx- Further, we denote 
the pulled-back metric X^gab of the physical spacetime A4x hy gab, as mentioned 
above. Though we have to keep in our mind that we are treating perturbations in 
the single gauge choice when we treat the derivative operator Va and the pulled-back 
physical metric gab on the background spacetime Mq, there is no confusion in the 
development of the perturbation theory if we treat perturbations only in the single 
gauge choice Xx- 

Since the derivative operator Va (= X*\/a {X~^) ) may be regarded as a deriva- 
tive operator on the background spacetime that satisfies the property (|2-30p . there 
exists a tensor field C^^j^ on the background spacetime A^o such that 

VaUJb = ^aU^b - C'ab^c, (2-31) 

where cUc is an arbitrary one- form on the background spacetime Mq. Prom the 
property (j2-30p of the covariant derivative operator Va on Ai, the tensor field C^^f^ 
on A^o is given by 



Cab = ^9 {^aQdb + ^b9da " dQab) , 



(2-32) 
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where g""^ is the inverse metric of Qah, i-e., gacO'^^ = (^a^- We note that the gauge 
dependence of the derivative Va as a derivative operator on M.o is included only 
in this tensor field C"^^^. The Riemann curvature R^bc'^ physical spacetime 

Aix, which is pulled back to the background spacetime A^o, is given by the Riemann 
curvature R^f^^ on the background spacetime M.q and the tensor field C^^^^ as follows: 

^abc = ^abc ~ 2V[^C"^^]^ + 2C^^[„C"^j,]g. (2-33) 

The perturbative expression for the curvatures are obtained from the perturbative 
expansion of Eq. (j2-33p through the perturbative expansion of the tensor C^^^ defined 
by Eq. (|M2]l . 

The first- and the-second order perturbations of the Riemann, the Ricci, the 
scalar, the Weyl curvatures, and the Einstein tensors on the general background 
spacetime are summarized in KN2005. We also derived the perturbative form of 
the divergence of an arbitrary tensor field of second rank to check the perturba- 
tive Bianchi identities in KN2005. In this paper, we only present the perturbative 
expression for the Einstein tensor. 

We expand the Einstein tensor on the physical spacetime M.\ as 

G," = G," + A(i)G,^ + ^A^Wg,^ + 0(A3). (2-34) 

As shown in KN2005, each order perturbation of the Einstein tensor is given by 

WG,^ = «^,''[?^] + £xG„^ (2-35) 

^'^G,* = ^^'^Ga + ^^^Ga [n,n] + 2£x^''^Gj' + {£y - £jc}Gj', (2-36) 

where 

(1)^." [A] := ^')eJ> [A] - \5j'^'^S,'^ \A] , (2-37) 

[A, B] := (2)i7,^ [A, B] - ^J.^^^^r/ [A, B] , (2-38) 

[A] := -2V^^H/'^ [A] - A^^R^c, (2-39) 

[A, B] := 2RadBj'A^^^ + 2H^/^ [A] [B] + 2i7[/^ [B] [A] 
+2A/Vy,H^^^- [B\ + 2i?/V[,F,j''- [A] 

+2Aj^Vy,H^^ [B] + 25, -'^ [A] , (2-40) 



and 



H^,^[A]:=V^,A,f--V^Aab, (2-41) 



1, 

2 

Habc [A] := gcdH^b " [A] , HJ'^ [A] := g'"'H^/ [A] , 
H,'AA]:=gcdH,'''[A]. (2-42) 



We note that (^^^^^ [*] and (2)^^^ [*, ^] in Eqs. (IM5]> and (fM6]l are the gauge in- 
variant parts of the perturbative Einstein tensors, and Eqs. (j2-35p and (j2-36p have 
the same forms as Eqs. (|2-26p and (j2-29p . respectively. 
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We also note that (^^^a^ [*] and (2)^^'' [*, *] defined by Eqs. (l2^ - (f24n] ) satisfy 
the identities 

Vji)^,« [A] = -H.^'^iA] G,^ + H,^^[A] (2-43) 

'i/fead [5] A'^'^ + Hbad [A] S^^) G,'' 
+ fi/,„, [5] ^'^'^ + H,ad [A] B'^'') (2-44) 



for arbitrary tensor fields Aab and Bab, respectively. We can directly confirm these 
identities, and these identities guarantee the first-order and second-order perturba- 
tions of the Bianchi identity VbGa — 0; respectively, as shown in KN2005. These 
identities are also useful when we check whether the derived components of ^^^Qa 



* 



and are correct. 

Finally, we consider perturbations of the Einstein equation of first and second 
order. First, we expand the energy-momentum tensor as 

= T/ + A(i>r/ + ^A2(2>r/ + 0(A3). (2-45) 

Following the definitions ()2-26p and (j2-27p of gauge invariant variables, the gauge 
invariant variables " and '^^•'T^ " for the perturbations '•^^a'' ^^^a'' °f 
energy-momentum tensor are defined by 

{lyj^ a ._ {1>^ a _ ^^rj.^ (2-46) 
(2)^^a _{2>^a_2£^(l>^a_|^^_^2^| j.^a_ (2-47) 

In ^4.31 we show that these definitions of the gauge invariant part of the first- and 
second-order perturbation of the energy-momentum tensor are appropriate in the 
cases of both a perfect fluid and a single scalar field. Further, we impose the per- 
turbed Einstein equation of each order, 

(i)G„'' = 8^G Wt/, = 87,(- (2)j.^fe^ (2-48) 

Then, the perturbative Einstein equation is given by 

(1)^/ [7i] = SvtGWt/ (2-49) 

at linear order and 

(1)0,^ [£] + (2)^^^ [n, n] = SttG (2)t, * (2-50) 

at second order. These explicitly show that, order by order, the Einstein equations 
are necessarily given in terms of gauge invariant variables only. Therefore, we do 
not have to consider the gauge degree of freedom, at least at the level where we 
concentrate only on the perturbed Einstein equations. 

We have reviewed the general outline of the second-order gauge invariant per- 
turbation theory. Within this general framework, we develop a second-order cosmo- 
logical perturbation theory in terms of the gauge invariant variables. 
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§3. Cosmological background spacetime 

Here, we consider ttie background spacetime for cosmological perturbation the- 
ory. The background spacetime considered here is a homogeneous, isotropic universe 
that is foliated by the three-dimensional hypersurface which is parameterised 

by 77. Each hypersurface of S{ff) is a maximally symmetric three-space, and the 
spacetime metric of this universe is given by 

gab = a^iv) {-id'n)aidr])b + 1ij{dx')a{dx^)b) , (3-1) 

where a = a{r]) is the scale factor, jij is the metric on the maximally symmetric 3- 
space with curvature constant K, and the indices i,j, k, ... for the spatial components 
run from 1 to 3. Depending on the behavior of the scale factor a, this metric (j3-ip 
can represent a Friedmann-Robertson- Walker universe or a de Sitter spacetime. In 
terms of the coordinate system in which the spacetime metric is given by (|3-ip . the 
components of the Christoffel symbol of this background spacetime are given by 

r^^ = Ti, -Tj^ = = r^^, r^j = H-jij, ri^ = T-c^fJ ^ r^j = ^ ^r^j, (3-2) 

where 7i := 9^a/a, ^^^Fjj is the Christoffel symbol associated with the three metric 

7jj, and 7j"' = JikJ^'' is the three-dimensional Kronecker delta. These components 
of the Christoffel symbol are useful when we write down the components of the four- 
dimensional covariant derivative of tensors in terms of the derivative with respect to 
rj and the three-dimensional covariant derivative associate with the metric 'jij. 

Since "jij is the metric on the maximally symmetric 3-space with the curvature 
constant K, the curvatures associated with the metric "jij are given by 

^^%ki = 2K^k[ilj]i, ^^% = 2Kji,, (% = 6i^. (3-3) 

These are useful when we calculate the components of the perturbative curvatures 
in terms of the 3-|-l decomposition, as in the metric (jS-ip . The four-dimensional 
background curvature tensors are also necessary to calculate the components of the 
perturbative curvatures. These are given by 

Rab = -3dnn {dri)^ (dr^), + {d^^H + 2n^ + 2K) 7,^, (3-4) 



(3-5) 



R," = ^d,n [dr,)^ + 1 + 2^2 + 2K) 7,^ 

R:=R^'' = ^{d^n + n'^ + K] , (3-6) 



Ga '■— Ra - 2^a^R 

b 



= -^[^' + K] idr^)a {^) - ^ [2dr,n + n' + K] 7,^ (3-7) 
where jab = lij{dx'-)a{dx^)b and 7^'' = j-\dx'-)a{d/dx^)''. 
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To study the Einstein equation for this background spacetime, we introduce the 
energy-momentum tensor for a perfect fluid, 

= euau'+p{6,'' + uau'') (3-8) 
= -e(dry), +p^J, (3-9) 



b 



where we have used 

We also consider the energy-momentum tensor for the scalar field, which is given by 

= Va^VV - (Vc^VV + 2V{ip)) (3-11) 



+ (y^{d,^f-y{^))la\ (3-12) 



where we have assumed that the scalar field (p is homogeneous, i.e., (p = (p{rj). Com- 
paring (|3-12p with (j3-9p . the energy density and the pressure for the homogeneous 
scalar field are given by 

'^ = ^idvV>? + V{ip), p = ^id,pf-V{^). (3-13) 

The Einstein equations Gj' = 8itGT^ ^ for this background spacetime filled with 
a perfect fluid are given by 

7i'^ + K = —-a^e, 2d^n + 71^ + K = -SiiGa^p. (3-14) 

In the derivation of the perturbative Einstein equations, the equation 

Ti^ + K -d^rL = ^^Ga^{e+p) (3-15) 

is also useful. Of course, there is an equation for the energy conservation of the 
matter fields, and this equation gives the behavior of the energy density in the scale 
factor if we apply an appropriate equation of state for the matter field. This equation 
is consistent with the two equations in (j3T4p . 

Further, in the case of the single scalar field model, the Einstein equations are 
given by 

n' + K = + y{^)) , (3-16) 
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through the relations (I3-I3p . We also note that the equations p-16p and (j3-17p lead 
to 



rC^ + K - driH = AiTGidriipf . (3-18) 

Equation (j3-18p is also useful when we derive the perturbative Einstein equations. 
Equations (j3-16p and (j3-17p are often used to investigate the inflationary scenario. 
Actually, in the situation that the potential term of the scalar field is sufficiently 
larger than its kinetic term, the spacetime is approximately a de Sitter spacetime, 
and this situation may be realized in the very early universe. Hence, the background 
spacetime described by Eqs. (|3-16p and (j3-17p also includes inflationary universes, 
and the second-order perturbation theory developed below is also applicable to the 
inflationary universe. 



§4. Gauge invariant variables of cosmological perturbations 



Now, we develop the second-order perturbation theory with the cosmological 
background spacetime in ^ within the general framework of the gauge invariant 
perturbation theory reviewed in ^ The important step when we apply the above 
general framework of the gauge invariant perturbation theory is to confirm that the 
assumption for the decomposition ()2-22p of the linear-order metric perturbation is 
correct. This confirmation is accomplished in ^4.11 Hence, the general framework 
reviewed in ^ is applicable. Applying this framework, we define the second-order 
gauge invariant variable of the metric perturbation in N4.2l and of the matter pertur- 
bations in 



4.1. First-order metric perturbations 

On the background spacetime discussed in fj3l we consider the metric perturba- 
tion to be that given in Eq. (j2-2ip . To show that the assumption for the decompo- 
sition (j2-22p of the linear-order metric perturbation is correct, we first consider the 
components of the linear-order metric perturbation in the coordinate system (|3Tp . 



Kb = Krjidr])a{dr])b + 2/i^j((ir/)(a((ix*)fc) + hij{dx')a{dx^)b- (4-1) 



Because components belonging to different groups are coupled through contraction 
with the metric tensor and the covariant derivatives in the Einstein equations, the 
grouping {hrirj-, h^/i, hij} is not so useful. Instead, using the fact that the three mani- 

''Til 



fold {I!{r]), a^jab) is maximally symmetric, we further decompose the vector /i^j and 



tensor hij as 



hrji = Dih(yL) + h(y)i, D'h(y)i = 0, (4-2) 

hij = a^h(L)-iij + a^h(T)ij, h^T)\ ■= l^^hT)ij = 0' (4-3) 

h{T)ij = (^DiDj - ^lij^^ hiTL) + '^D(^ih(^TV)j) + h{TT)ij, (4-4) 

D'h^TV)i = 0, D'h(^TT)ij = 0, (4-5) 
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where Di is the covariant derivative associated with the metric jij, and A = D^Di. 
It is weh-known that the linearized Einstein equations on a homogeneous isotropic 
universe, whose metric is given by Eq. (j3-ip . can be decomposed into groups that 
each contains only variables belonging to one of the three sets {/i,jr;, h(VL)i ^{L)-> ^tl}-, 
/i(2-y)^}, and /i{TT)ij Variables belonging to these sets are called scalar- 
type, vector-type and tensor-type variables, respectively. This segregation of vari- 
ables is due to the fact that the metric tensor a^^ab is the only non-trivial tensor on 
the maximally symmetric space, and as a consequence, the tensorial operations on 
hab to construct the linearized Einstein tensors preserve this decomposition. 

In the linear perturbation theory, the covariant derivatives are always combined 
into the Laplacian A in the linearized Einstein equations after the decompositions 
(l4-2[ )- (j4-4p . because the metric tensor a^^ab is the only non-trivial tensor on the 
maximally symmetric space. Thus, the harmonic expansion of the perturbation 
variables with respect to the Laplacian is also useful in the linear perturbation theory. 
However, in the second-order perturbation theory, mode-mode coupling occurs due 
to the non-linearity of the Einstein equations. For this reason, we do not apply the 
harmonic expansion of the perturbation variables with respect to the Laplacian in 
this paper, though the harmonic expansion should be also useful after the nonlinear 
terms in the second-order Einstein equations are clarified. Instead, we assume the 
existence of some Green functions, as explained below. 

To clarify the uniqueness of the decompositions (j4-2|) - (j4-4|) . we consider the 
inverse relations of Eqs. (j4-2p - (j4-4p . To do this, we first note that the commutator 
of the divergence and the Laplacian is given by 



D'Ati - AD% = D' i^^^^Rijt^j = 2KDHi, (4-6) 
D'AU, - AD%, = D' (^^%He, + + ^^%kjiD^f' 



AK[D\,^^-\D,t^] (4-7) 



for any tensor fields ti and tjj, where we have used Eqs. (j3-3p . These relations are 
also useful when we write down the components of the perturbative Einstein tensor. 
They show that on the maximally symmetric space, the Laplacian preserves the 
transverse condition 

D^Ah^v)j = + 2i^) D^h^v)j = 0. (4-8) 
The inverse relations of the decompositions ()4-2p - (j4-4p are given by 

h^VL) = ^'^D^hrjj, (4-9) 
h^V)^ = Ki-D,A-^D^h^,, (4-10) 

hL) = ^K\ (4-11) 
hT)ij = ^ (hij - ^/i^T^n , (4-12) 
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Vl) = \{^ + 3K)-'A-'D'D^h(^r)ij, (4-13) 
h^TV)^ = + 2Kr^D%T)^k " + 2K)-^ D,A-^ D' h(^T)kU (4-14) 

-2D(,(Z\ + 2A')-^D'=V)j)fc 

+2D(,(Z\ + 2K)-'D^^A-'D''D'h^T)ki- (4-15) 

Equations (j4-9p - (|4-15p show that there should be exist Green functions of the op- 
erators A, A + 2K, and A + 3K to guarantee the one to one correspondence of 
the decompositions of the set {hr/r^, hrji, hij} and the sets { {/i^r,, /i(y^), /i^x,), /itl}, 
., /i(Ty)^}, /i(Tr)ij }• Actually, these Green functions exist if we specify the do- 
main of the perturbations, for example L^-space on S{r]) with appropriate boundary 
condition. Therefore, we assume the existence of these Green functions in this paper. 
By this assumption, any tensor that belongs to the kernel of any of the operators A, 
A + 2K, and A + 3K is excluded from consideration. For example, a Killing field 
Vi on the three-dimensional hypersurface 2J{ri), which satisfies the Killing equation 
^(i'^j) ~ ^' belongs to the kernel of the operator A + 2K, since we can easily confirm 
(A + 2K)vi = from the Killing equation. If it is necessary to investigate such 
tensors as the perturbative mode, separate treatments are necessary. Because the 
treatment of these exceptional modes is beyond the scope of this paper, we ignore 
all modes that belong to the kernels of the operators Z\, Z\ + 2K, and A + S-fT. 

Now, we consider the decomposition (j2-22p of the linear-order metric pertur- 
bation and show that the decomposition (j2-22p is valid in the case of cosmological 
perturbations, though this decomposition is merely an assumption in the general 
framework reviewed in ^2.2[ 

To accomplish the decomposition (j2-22p . we consider the gauge transformation 
rule (j2-19p . which is given by 

jAab - xKb = £^gab = 2V(a^b) (4-16) 

for linear-order metric perturbations. In Eq. (j4-16p . the generator of the gauge 
transformation is an arbitrary vector field on the background spacetime A4q. We 
decompose the generator in terms of the 3+1 decomposition as 

= Udv)a + Udx%, (4-17) 
and, further, the component as 

C^ = D.i^L) + C(T)i, D'i[T)i = 0. (4-18) 

In terms of the 3+1 decomposition, the gauge transformation rules (j4-16p are given 
by 

:ihrrn- xhr,ri = 2{dr^-'H)ir,, (4-19) 
A* - xhr1^ = + {d^ - 2n) Ci, (4-20) 
^ij - xhi, = 2L»(,e,) - 2n-(i,i^,. (4-21) 
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Furthermore, following the decomposition (j4-18p of the component ^j, the gauge 
transformation rules (j4-20p and (j4-2ip are obtained in terms of ^rj, ^(l), and C{T)i as 

yhr,i - xhrji = A { {drj - 2n) + C„} + (5r, - 2H) i^T)^, (4-22) 
■^ij - xhij = 2 Qz\C(i) - ■jij 

+2 (^D,D, - ^7,,Z\^ + 2D(^iC\:^K (4-23) 

Since the tensorial decomposition (j4-2p - ()4-4p has the inverse relations ()4-9p - ()4-15p . 
the gauge transformation rules (j4-22p and (j4-23p lead to those for the metric pertur- 
bations /i(yL), /i(y)i, /i(rL), /i(ry)i, and /i(TT)ii as follows: 

^{VL) - xh(yL) =iv + - 2W) kL) ' (4-24) 
:^{v)i - xh{v)i = (dr^ - 271) S,(^T)i, (4-25) 

aV(L) - a^xhiD = -2^., + ^/^^(L), (4-26) 
a^:^(TL) - a^xh{TL) = 2i[L), (4-27) 
a^:^{TV)i - d^xh{TV)i = i{T)i-, (4-28) 



a 



-'^yh(^TT)ij - a^xhiTT)ij = 0. (4-29) 



Inspecting the gauge transformation rules (|4-24p ~ (|4-29p . together with Eq. (|4-19p . 
we find gauge invariant and variant variables of the first-order metric perturbation. 
We first construct gauge invariant variables. 

First, equation ()4-29p shows that the transverse traceless part /i(TT)jj is itself 
gauge invariant. We denote this transverse traceless part by 

Xij ■= l^{TT)ij, Xij = Xji, X i=0, D Xij = 0. (4-30) 

(1) 

The transverse traceless tensor Xij has two independent components and is called 
the "tensor mode" in the context of cosmological perturbations. This transverse 
traceless part of the metric perturbations is well known as gravitational waves on a 
homogeneous isotropic universe. 

Next, we consider to the gauge transformation rules (j4-25p and (j4-28p . From 
these gauge transformation rules, we can easily see that the variable defined by 

% := h(v)i - (dri - 271) [a^\TV)i) 

= h{v)i - a^d^h(TV)i (4-31) 

is gauge invariant. The gauge invariant variable z/j is called a "vector mode" in the 
context of cosmological perturbations. It satisfies the equation 

% = (4-32) 
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from the divergenceless property of the variables and h(^xv)i- Equation (I4-32P 

(1) 

imphes that the vector mode includes two independent components. 

In addition to the vector and tensor mode of the perturbation, there are two 
scalar modes in the linear-order metric perturbations hat- To see this, we first 
consider the gauge transformation rules (l4-24p and (I4-27P . From these transformation 
rules, the variable defined by 

■■= h{vL) - ^ (5»7 - 2?^) (a^/i(rL)) = h^vL) - ^a'^dr,h(^TL) (4-33) 

is transformed as 

yXr^ - xXr, = yh(^vL) - xh(vL) - -^{dr)- 2H) [a^ {yh{TL) - xh{TL))) 

= dr,i(^L) + - 27ie(L) -\^dr,- 271) (2e(i)) . 

= e,- (4-34) 

Using and inspecting the gauge transformation rule (j4T9p . we easily find that 
(1) 

the variable ^ defined by 

- 2a^ ''S := hrjr^ - 2 {d^ - 7i) Xr, (4-35) 

is gauge invariant. Further, from gauge transformation rules (j4-26p . (j4-27p . and 

(1) 

(|4-34p . the variable ^ defined by 

- 2a^ ^^■.= {h^L) - ^Ah^^L)^ + 2WX^ (4-36) 

(1) (1) 

is gauge invariant. The two scalar functions and ^ are called "scalar perturba- 
tions" in the context of cosmological perturbations. 

Thus, we have six components of gauge invariant variables: two components of 
the tensor mode, two components of the vector mode, and two scalar modes. Since 
the metric perturbation hab has ten components, there are four remaining compo- 
nents, which are the components of the gauge variant part of the metric perturbation. 

Because we already have all gauge invariant variables, we can specify the variant 
part Xa of the metric perturbation hab as in Eq. (j2-22p . Using the gauge invariant 
variables ^, i^i, and Xij, the components of the metric perturbation hab are given 
by 

hrjrj = -2a^ ''S +2 {dr, - 7i) Xr,, (4-37) 
hjji = +a'^dr,h(^TV)i + Dih^vL)-, (4-38) 

hij = -2a 'If jij + a Xij +a DiDjh^TL) - 27iXnJij + 2a D^ih(j^v)j)- (4-39) 
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On the other hand, in terms of the 3+1 decomposition, the components of (I2-22P are 
given by 

hr,,,=n^r, + 2{dr,-n)X^, (4-40) 

h^, = n^, + DiX^ + d^Xi - 2nXi, (4-4i) 

h,j = Hij + 2Z?(,Xj) - 2n-f^jX^. (4-42) 
Here, Tiab is gauge invariant, and its components should be identified by 

2 2 (1) 2 2 (1) /X 

Hr,n ■= -2a T^^j := a i^i , Hij := -2a ^ jij + a Xij • (4-43) 
Then, we obtain the equations for X^. 

2{dr,-n)x^ = 2{dr,-n)x^, (4-44) 

D,X„ + {d^ - 2H) Xi = a^^^h^TV)^ + Dih^VL) , (4-45) 

2D^iXfj - 2H^ijX^ = a^DiDjh(^rL) - ^HX^j^^ + 2a^D(^,h(^TV)j)- (4-46) 

Equation (|4-44p yields 

Xrj = Xrj + aC^, (4-47) 

where (7^ is a scalar function satisfying the equation drjCrj = 0. Substituting (14-470 
and (|T33l) into ([¥I5|) . we obtain 

= a2 + ^ A/i(TL)) - DiC^a^ a^Q, (4-48) 

where Cj is the vector field satisfying the condition = 0. Substituting (|4-47p 
and (I4-48P into (|4-46p . we obtain 

a^D^iCj) - DiDjCr,a^ j ~^ ~ '^^^j"^'? = °- ^^'^^^ 

Thus, we have found that the gauge variant part of the metric perturbation X^ is 
given by 

Xr, = X^ + Cr, = \vL) - \a^drh(TL) + Cr,, , (4-50) 

Xi = (h^TV)i + ^A/i(TL)) + Ci, (4-51) 

where and Cj are defined by 

:=aCr^, (4-52) 

Ci := -D^Cy / ^ + a^Ci, (4-53) 



a 



and drjCri = = S^Cj. 
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Through Eqs. (j4-52p and (l4-53p . with the constraint (j4-49p . it is easy to confirm 
that the vector field defined by 

Ca:=Cr,{dr^)a + Ci{dx')a (4-54) 

is a Killing vector on the background spacetime M.q. Actually, it is readily shown 
that 

2V(^C7^) = 2{d^- H) Cr, = 0, (4-55) 

due to the definition (j4-52p and dr^Cn = 0. Further, from the definition ()4-53p . we 
can easily confirm that 

2V(^Ci) = ^^C^ + AC^ - 2rLCi = 0. (4-56) 
Finally, the constraint (14-49P leads to 

2V(,C,) = 2D(,C,) - 2rL-iijCr, = 0. (4-57) 

Thus, we have 

V(„Cb) = 0; (4-58) 

i.e., Ca defined by Eq. ()4-54p is a Killing vector. Hence, we have been able to specify 
the gauge variant part Xa of the linear-order metric perturbation as 

Xa:=Xr,{dri)a+Xi{dx')a. (4-59) 

The relation between the components of the gauge variant part, Xa, and the com- 
ponents of the linear-order metric perturbation, h(y^, ^(tl); h(xv)i, is determined 
up to the degree of freedom of the Killing vector field. Since Xa contributes to the 
metric perturbation as in Eq. (l2-22p . the Killing vector field Ca in Eqs. ()4-50p and 
()4-5ip does not contribute to the metric perturbation. 

Further, because we do not consider the kernels of the operators Z\, A + and 
A + 3-fr as the domain of the perturbations, the gauge variant part Xa of the metric 
perturbation does not have the degree of freedom of the Killing vector field; i.e., the 
gauge variant part Xa is determined without ambiguity. To satisfy Eq. ()4-49p for 
any scale factor a, the components and Ci should satisfy the equations 

Cr, = 0, ^(,(7,) = 0. (4-60) 

This implies that Cj is a Killing vector on a three-dimensional hypersurface S{ri). 
As commented just after Eqs. (|4-9p - (|4-15p . it is easily shown that {A + 2K)Ci = 0. 
Therefore, the vector field Ci does not belong to the domain of perturbations con- 
sidered here. Thus, we have Cj = 0. Hence, we conclude that we can determine the 
gauge variant part Xa without ambiguity. Of course, it might be possible to include 
the Killing fields in our consideration by extending the domain of the perturbations. 
Separate treatments are necessary to do this, as mentioned above. 

Finally, we check the transformation rules for the vector field Xa under the 
gauge transformation X\ — > y\. Because the component Xr^ of the vector Xa is 
transformed as 

yX^ - xXr, = Cv, (4-61) 
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as noted in Eq. (l4-34p . From Eq. (l4-5ip with Cj = and the gauge transformation 
rules ()4-27p and (j4-28p . the gauge transformation rule of the component Xi defined 
by Eq. (j4-5ip is given by 

2 2 / 2 2 ^ 

yXi - xXi = a yh(TV)i - a xh(TV)i + [a Ji{tl) - a xh(TL)) 

= ^iT)^ + Di^iD = ^i- (4-62) 

Together with the transformation rules ()4-6ip and (j4-62p . the vector field Xa defined 
by Eq. (|4-59p is transformed as 

yXa - xXa = (a- (4-63) 

This shows that Xa is the gauge variant part of the metric perturbation in the 
decomposition ()2-22p . 

Thus, we know the procedure to find the gauge invariant variable Tiab and the 
gauge variant variable Xa in the case of cosmological perturbations. We have con- 
firmed the important premise of the general framework of the second-order pertur- 
bation theory reviewed in ^ Hence, we can apply this general framework for the 
second-order gauge invariant perturbation theory presented in KN2003 and KN2005 
to cosmological perturbations. 

4.2. Second order metric perturbations 

Here we consider second-order metric perturbations. We expand the metric gab 
on the physical spacetime Mx as Eq. (j2-2ip . According to the gauge transformation 
rule (j2-20p . the second-order metric perturbation lab is transformed as 

^iab - ?^ab = '^£((i^xhab + |'^5(2) + "^lij} 5a6 (4-64) 

under the gauge transformation <Px = {Xx)~^ ° y\ '■ ^^x ^ 3^a- As shown in ^4.H 
the first-order metric perturbation hat is decomposed in the form (j2-22p . Using this 
important fact, as shown in KN2003, the second-order metric perturbation lab can 
be decomposed as Eq. ()2-24p . Here, we demonstrate this. 

Inspecting the gauge transformation rule (|4-64p . we first introduce the variable 
Lab defined by 

Lab ■= lab - '^£xKb + £x9ab- (4-65) 

Under the gauge transformation <Px = (-^a)^^ ° '■ '^x ^ Vx, the variable Lab is 
transformed as 

yLab — xLab = J^afe " '2£yX:^ab + £yX9ab 

-?S'ab + '^£xXxhab - £\x9ab 
= 2£^^^^xKb + {^£(^2) + £\r^9ab 

-2£^X+5(i) {xKb + £^^^^9ab^ + £l.x+i^i^9ab 

-^iab + '2£xXxhab - £\x9ab 
— £a9abi 



(4-66) 



(4-67) 
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where 

<T'':=it2) + [ki).Xr- (4-68) 

The gauge transformation rule (|4-67|) is identical to that for a linear metric perturba- 
tion. Therefore, we may apply the procedure to find the gauge invariant and variant 
variables of linear-order metric perturbations to the decomposition of the compo- 
nents of the variable Lab- Arguments completely analogous to those used in the case 
of linear-order metric perturbation show that the variable L^b can be decomposed 
as 



Lab = Cab + £Y9ab, (4-69) 

where Cab is the gauge invariant part of the variable 1/^^, or equivalently, of the 
second-order metric perturbation lab, and Y"- is the gauge variant part of Lab, i-e-, 
the gauge variant part of lab- Under the gauge transformation <P\ = [X\)^^ °y\-, the 
variables Cab and are transformed as 



yCab - xCab = 0, y^a - :^a = CTa, (4-70) 

respectively. Thus, we have reached the decomposition (|2-24p of the second-order 
metric perturbation lab into the gauge variant and gauge invariant parts. Following 
to the same argument as in the linear case, the components of the gauge invariant 
variables Cab are given by 

Cab = -2a^ y {dr])a{dr])b + 2a^ % {dri)(^a{dx')b) 

W (^-2 -fij+ %j ^ {dx')a{dx^)b, (4-71) 

(2) (2) 

where v i and Xij satisfy the equations 

■ (2) • (2) 

D'u^=^, X\ = 0, D'x^J = ^- (4-72) 

, (2) (2) . , 

The gauge invariant variables <p and are the scalar mode perturbations of second 

(2) (2) 

order, and Vi and Xij are the second-order vector and tensor modes of the metric 
perturbations, respectively. 

4.3. Matter perturbations 

Since we have obtained the first- and the second-order gauge variant parts, Xa 
and Ya, of the metric perturbation, we can define the gauge invariant variables for 
an arbitrary field Q, except for the metric by following the definitions (|2-26p and 
(IJ27ll .^^) Here, we consider the first- and second-order gauge invariant variables for 
the perturbations of the perfect fluid components, the single scalar field, and their 
energy-momentum tensors. 
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4.3.1. Perfect fluid 

First, we consider the perturbation of a perfect fluid. As shown in Eq. ()3-8p . the 
total energy-momentum tenor of the fluid is characterized by the energy density e, 
the pressure p, and the four-velocity u": 

fJ = {e + p)uau'' + pSj. (4-73) 

Of course, this energy-momentum tensor is the representation on the physical space- 
time M\, but we can regard this equation to be the representation on the background 
spacetime A^o which is pulled back by an appropriate gauge choice Xx- The back- 
ground value of the energy momentum tensor (j4-73p is given by Eqs. (|3-8p - (|3T0p in 
^ In addition to the components of Eq. (j4-73p . we may also include the anisotropic 
stress as non-diagonal space-space components of the energy-momentum tensor, as 
phenomenology. Because we can always extend our arguments to those including 
the anisotropic stress, we ignore it here, for simplicity. 

Now, we consider perturbations of the fluid components of the energy-momentum 
tensor (f¥73l) : 

e:=e + X^l\^X^^f +0{X% (4-74) 
p-p + AP+^A^ pVo(A3), (4-75) 

Ua ■.= Ua + X (Ua) +-X^ (n„) +0{X^). (4-76) 

The fluid four-velocities Ua on the physical spacetime and Ua on the background 
spacetime satisfy the normalization condition 

r^'UaUa = g^'^UaUb = -1. (4-77) 

The perturbative expansion of the normalization conditions (j4-77p gives the con- 
straints for the components of the first- and second-order perturbative four-velocities 

(1) (2) 

{ua) and {ua)- The perturbative expansion of the normalization condition (|4-77p to 
second order gives the normalization condition at each order: 

{Ua) = -h'^^aUb, (4-78) 
{Ua) = h'^'ua (Uh) -g'''' {Ua){ub) -h^'K UaUb + -l^^^UaUb. (4-79) 

We also consider the perturbation of the four- velocity as 

(1) 1 ^ (2) 

= u'' + X (u") +-A2 (n'^) . (4-80) 
The perturbative expansion of the equation leads to 

{u-) = g-' {ub) -h^'ub, (4-81) 

(2) (2) (1) 

(u") = g"'' (ub) -2h''^ (ub) +(2/i'^^/i,* - r'')ub. (4-82) 
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Further, the first-order perturbation (I4-78P of the normalization condition (j4-77p is 
given by 

(1) (1) 

ua + n'^ (na)=0. (4-83) 

Next, we define the gauge invariant variable for the perturbation of the fluid 
components e, p, and Ua- Following the definitions (|2-26p and (j2-27p of the gauge 
invariant variable for an arbitrary matter field, we define the variables 

£ := V -£xe, (4-84) 
(1) (1) 

V:= P -£xP, (4-85) 

(1) (1) 

:= (Ua) -£xUa, (4-86) 

f ■.= ^e -2£x^e -{£y - £x}^, (4-87) 
V := P -2£x P -{£y-£x}p, (4-88) 

(2) (2) (1) ^ 

Ua:={Ua)-2£xUa-{£Y-£x}ua, (4-89) 

where the vector fields Xa and Ya are the gauge variant parts of the first- and second- 
order metric perturbations, respectively, which were defined in ^S ]4.2l and 14. li 

The first-order perturbation (I4-78P of the normalization condition (j4-77p is given 

by 

S = lUabu'^u'' - £x Q«"«a) (4-90) 
= \nabu''u'', (4-91) 
while the second-order perturbation (j4-79p of Eq. (j4-77p is given by 

(2) (1) (1)(1) 1 

Ua = 2nabU''g'"' -5"' UaUb -nacHdbg''''UaUb + -CabU'^u'' 



-2£x iu'' {ua) -iKbu'^v!' 



'£y - £l) (^«an«) (4-92) 



(1) 1 

= 2nabu''g^' -5"' UaUb -Uacndbg'^'UaUb + -llabu''u\ (4-93) 

where we have used Eqs. ([¥77]) and ([¥78]) . We note that (|¥90l) and (|¥92]) have 
the same forms as definitions (|2-26p and (j2-27p of the first- and second- order gauge 
invariant variables for an arbitrary tensor field, respectively. These are natural re- 
sults, because Eqs. ()4-90p and (j4-92p are the results of the first- and second-order 
perturbative expansions of the variable (l/2)n"Ma. 

We also decompose the first- and second-order perturbations of the four-velocity 
u", which are given by (j4-8ip and (|4-82p . respectively, into gauge invariant and 
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variant parts: 

(n'^) = g""^ Ub -n^^Ub + (4-94) 
. (2) (1) 

+2£x 1^5"' (Si'b) -/i'^'n,^ + (i:y - £l) (u-) . (4-95) 

We note that these expressions have the same forms as Eqs. (j2-28p and ()2-29p . 

Next, we consider the expansion of the energy-momentum tensor ()4-73p . Substi- 
tuting the expansion (l4-74p - ()4-76|) of the fluid components e, p, and Ua into (14-730 . 
we have the perturbative form of the energy-momentum tensor: 

f/=:r/ + A«(r/) +iA2{2)(r/) +0(A3), (4-96) 



where 



{Ta ') = (^^^ + uau' + {e + p) Ua [u") + {e + p) {Ua) + P K\ (4-97) 

, , /(o) (2)\ (2) (2) 

(2) (r, " j = r + p J uav" + {e + p) Ua {u') +{e + p) {Ua) u'+ V dj' 



(1) (1) 

+2(e + p) K)(n^) . (4-98) 

Further, we consider the decomposition of the perturbed energy-momentum tensor 
into its gauge invariant and variant parts. Using the definitions (|4-84p ~ (j4-86p of 
the gauge invariant variables of the first-order perturbations and Eq. (j4-94p . the 
first-order perturbation of the energy-momentum tensor of a perfect fluid is given by 



«(r/) =-A^yrJ + £xTa\ (4-99) 
where the gauge invariant part of the energy-momentum tensor is given by 

(i)t/:= (s+ V\uau''+ V Sa' 

( (1) (1) \ 

+ {e+p)[uaU^ -n^%cUa+ UaU^y (4-100) 

Similarly, using the definitions (|4-84p - (j4-89p of the gauge invariant variables of the 
first- and second-order perturbations and Eqs. (j4-94p and (j4-95p . the second-order 
perturbation of the energy-momentum tensor of a perfect fluid is given by 

(2) (T^b\ (2)^^ . + 2i:x« (r/) +{£y-£l]Tj, (4-101) 
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where 



(2>r^ b^o(^'JJ^\ (lib 

^ :=(£•+ 7? I UaU +2\£+'P\Ua\U -H u. 



/ (2) (1) ^ 



,'(1) (1) , W f , (1) 



+ {e + p)UaU^+ V (4-102) 

Here, again, we have seen that the perturbative expressions (I4-99P and ()4-10ip of 
the energy-momentum tensor have the same forms as Eqs. (I2-28P and ()2-29p . as 
expected. We also note that in the derivation of the expressions (|4-99p - (|4-102p . we 
did not exphcitly use any background values of the fluid component and the metric. 
Therefore, the expressions (j4-99p - ()4-102p are valid for any background spacetime. 
This implies that the definitions (j2-46p and (j2-47p of the gauge invariant variables 
of the perturbed energy-momentum tensor in §2.41 are appropriate in the case of a 
perfect fluid. 

Finally, we consider the perturbation of the equation of state for a fluid. In the 
generic case, the equation of state of a fluid is given by 

p = p{e,S), (4-103) 

which gives the relation between the pressure p, the energy density e, and the entropy 
S. In addition to the perturbative expansions (14-741) and (I4-75P for the energy density 
and pressure, we consider the perturbative expansion of the entropy: 

5 = S + A 5 +^A2 +0{X^). (4-104) 
Hence, the generic equation of state (j4-103p is expanded as 

,^x'i> + 1 ? = i, + A 'J' +iA' 5 + A ^' +iA^ 'I) (4.105) 

/ m . fdp (1) dp (1) 

+ '|2(.,SH!^''g(.,S)V (4.106) 

Thus, we obtain the equation of state of the first- and second-order perturbation of 
the fluid components: 
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(2) (2) 5p , ^, (ifd^p^ ^(i)(i) d'^p . ^. 

(2) dp i^f d^p 

+ S^{e,S)+S 0^(^,3). (4-108) 

In addition to the definitions of the gauge invariant variables ()4-84p - (j4-89p for 
the first- and second-order perturbations of the fluid components, we also define the 
gauge invariant variables for the entropy perturbations: 

(1) (1) 

S := S -£xS, (4-109) 

(2) (2) (1) 

S := S -2£x S -{£y- £\] S. (4-110) 

Substituting Eqs. (j4-84p . (j4-87p . and (j4-109p into the perturbations of the equation 
of state ()4-107p , we obtain the first-order perturbation of the equation of state of the 
fiuid in terms of the gauge invariant variables: 

V +£xP = H +£x€^ + ^ (^S +£xS 

dp W dp W ^ / , - , - N 

The right-hand side of Eq. (|4-llip has the same form as Eq. (|2-28p . as expected 
by considering the left-hand side of Eq. (j4-llip . Hence, we obtain the first-order 
perturbation of the equation of state in terms of the gauge invariant variables: 

(1) , (1) (1) 

V=ci £ +T S, (4-112) 

where we have 

2 9p dp 

'^'■=d-t ^-■=W ^'-'"'^ 
and Cs is interpreted as the sound velocity of the fiuid. The equation (|4-112p is the 
equation of state for the gauge invariant variables of the first-order perturbation of 
the fluid components. 

Next, we consider the second-order perturbation of the equation of state of the 
fluid in terms of gauge invariant variables: 

V +2£x P + {£y - £x) P 
= £^^{e,S)+£ Q^i^^S)+2£S^{e,S)+Sg^{e,S)+S ^(e,^) 
+2£x (^e ^(e, S)+ 5 ^{e, S)^ + £yp{e, S) - £lp{e, S). (4-114) 

The right-hand side of Eq. (j4-114p has the same form as expected by considering the 
left-hand side of Eq. (|4-114p . Then, we obtain the second-order perturbation of the 
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equation of state in terms of the gauge invariant variables: 

(2) (2) (2) (1)^ dT (1)^ 

V =cl £ +T S £ +2-^ £ S +^ S , (4-115) 

where we have used Eqs. (|4-113|) . 
4.3.2. Scalar field 

Here, we consider the energy-momentum tensor of the single scalar field (f, 

fj = V.c^vV - 1^,' (V,^V^^ + 2F(^)) , (4-116) 

where V{(p) is the potential of if. Because we are considering a perturbation theory 
on a homogeneous and isotropic universe, the scalar field must also be approximately 
homogeneous. Hence, the scalar field (p can be expanded as 

(^ = (^ + A(^i + ^A2(^2 + 0(A3), (4-117) 

where ip = '^{ri) is a homogeneous function on the homogeneous isotropic universe. 
The background field is the homogeneous part of the scalar field which drives 
the background homogeneous isotropic model, and \(pi\ -C \'f\ and |(^2| ^ I'^il are 
the first- and second-order perturbations of the scalar field (p, respectively. The 
energy-momentum tensor (|4-116p can also be decomposed into the background, the 
first-order perturbation, and the second-order perturbation as 

r/ = r/ + A « (r/) + (2) j^j.^^^ ^ ^(A^), (4-118) 

where (^^(T^'') is linear in the metric and the matter perturbations hab and (pi, and 
(2)^2T^6^ includes the second-order metric and matter perturbations, lab and (p2, and 
the quadratic terms of the first-order perturbations, (pi and hab- 

Expanding the metric as in Eq. (j2-2ip and the scalar field as in Eq. (j4-117p . 
the perturbations (^^(T^'') and ^'^\T^^) of the energy-momentum tensor (j4-116p are 
given by 

(t/) := Va^VVi - Va(^/i'"=V,^ + V,(^iVV 

-\^a (vcV^VVi - Vciph'^'^Vd'P + Vc(^iVV + ^'^il^) ' (4-119) 
(2) (j^J^ Va99VV2 - 2VaVP/i''^Ve(^i + Vaif (2h'''h/ - /^^) V^v? 

+2Va(^lV''(^l - 2Va^lh^^Vcip + Va(^2/'Vc(^ 



+2Vc(^iV^(^i - 2Vc(^i/i'^^Vrf(/5 + Vc^2VV 

+2(^2^ + 2((^i)2^ . (4-120) 
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According to the decompositions (j2-28p and (l2-29p . the perturbations of the 
scalar field ip at each order can be decomposed into the gauge invariant and variant 
parts as 

ipi =: ipi + £x^, (4-121) 
^2=:^2 + 2£x^i + {£y - £x)^, (4-122) 

where (pi and are the first- order and second-order gauge invariant perturbations 
of the scalar field. Through these gauge invariant variables, the perturbed energy- 
momentum tensor at each order can also be decomposed into the gauge invariant and 
variant parts. Substituting Eqs. (|4-12ip and (|2-22p into Eq. (|4-119p . the first-order 
perturbation (I4-119P of the scalar field is given by 

«(r/) =-A'>T,' + £xTj, (4-123) 

where 

is the gauge invariant part of the first-order perturbation of the energy-momentum 
tensor for the single scalar field. Through Eqs. (|4-12ip . (|4-122p . (p^ . and ([2^ . 
the second-order perturbation (j4-120p of the energy-momentum tensor is given by 

(2) (r/) =: (2)T, ' + 2£x« (t/) + {£y - £\) tJ, (4-125) 

where 



2' 

1 



+2Va(^iVVl - 2Vaipin^''Vcip + VaV92VV 



2 



+2ip2^ + 2ip{—^ \ . (4-126) 



The tensor '•^-'7^ * is the second-order gauge invariant part of the energy-momentum 
tensor for the single scalar field. Here again, we have seen that the perturbative ex- 
pressions (|4-123p and (j4-125p of the energy-momentum tensor have the same forms 
as Eqs. (|2-28p and (|2-29p . respectively, as expected. We also note that in the deriva- 
tion of the expressions (j4-123p - (|4-126p . we did not exiplicitly use any background 
values of the scalar field nor the metric. Therefore, the expressions (|4-123p - (j4-126p 
are valid for any background spacetime. This implies that the definitions (|2-46p and 
(j2-47p of the gauge invariant variables of the perturbed energy-momentum tensor in 
'> \2A\ are appropriate not only in the case of a perfect fluid but also in the case of a 
single scalar field. 
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§5. First-order Einstein equations 



In this section, we consider the perturbed Einstein equations of hnear order, 
(j2-49p . To derive the components of the gauge invariant part of the linearized Ein- 
stein tensor [Ti], which is defined by Eqs. (IM7D and (l2^ . we first derive the 
components of the tensor H^^^'^ \H], which is defined in Eq. (12-4ip with Aah = Ti-ab- 
Since the components of the gauge invariant part Tiab of the first-order metric per- 
turbation are given by Eq. (j4-43p . the components of the tensor H^^^[7{\ are as 
follows: 



H^^^ [H] = dr, 

H,^ " in] = A S +n % 



H., " [n] 



(5-1) 
(5-2) 



/(i) ^ (1) 1 (1) , , 

2nU + cp] +dr, ^\-fij- l/j) +- {dr, + 2H) X^3 , (5-3) 

(5-4) 



. (1) 

D' ^ + {dr, + n) v\ 



(1) 



1 



(1) 



(1) 



V,- + -dr, X/ 



(5-5) 



(1) 



(1) 



Hj,^ [H] = A 'i^ jkj - 2j\,Dj^ -Hjk, V* +Z)(, % . (5-6) 

The components of the tensors Habc [H], HJ"" [H], HJ"^ [H], and H""^" [H] are also 
useful when we derive the components of the gauge invariant parts ^^^QJ* and 
^"^^QJ* [Ti., Ti] of the perturbative Einstein tensor. These components are summarized 
in Appendix [Al 

Following to the definitions (j2-37p and (|2-39p of the gauge invariant part ^^^Ga^ [Ti] 
of the first-order perturbation of the Einstein tensor, its components are derived as 
follows: 



^'^Grj'in] 



1 
1 



-endr, + 2A + 6K) ^ -m^ } , 



2dr,D, ^ +2HA <P --{A + 2K) 



1 



(1) (1) 1 (1) 

-2 <^ 2dr,D' <F +2nD' ^ +- {-A + 2K + - Ad^Ti) 



(5-7) 
(5-8) 

(5-9) 



1 



AA - ^ 



+ 1 (-z\ + 2dl + mdr, - 2K) + {2ndr, + Adr,n + 2n^ + a) l^^j 7. 



1 



(1) 



(1) 



-^dr, i a' A V, 
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1 



(5-10) 



Straightforward calculations show that these components of the first-order gauge 
invariant perturbation ^^'^QJ' [H] of the Einstein tensor satisfies the first-order per- 
turbation (j2-43p of the Bianchi identity. This implies that we have derived the 
components (IPfll- ffFTOl) of [H] consistently. 

Together with the components of the gauge invariant part ^^^7^ ^ of the first- 
order perturbation of the energy-momentum tensor, the first-order Einstein equation 

()2-49p is given as equations for the gauge invariant variables ^ , "If , and x!j • We 
consider these equations for two cases, that in which the energy- momentum tensor 
is dominated by the single perfect fluid and that in which it is dominated by the 
single scalar field. 

5.1. Perfect fluid case 

Here, we consider the linearized Einstein equation of a homogeneous isotropic 
universe filled with a perfect fluid. 

(1) 

We first consider the components of the gauge invariant part lAa of the pertur- 
bative four-velocity of the fluid. Taking into account the perturbation (j4-9ip of the 

(1) 

normalization condition (l4-77p . the components of lAa are decomposed as 

(1) (1) / (1) . (1) 

Ua= -a <p {dr])a + A V + Vi {dx')a., V^= 0, (5-11) 

(1) 

where the T/-component of Ua are determined by Eq. (|4-9ip . Here, we note that 

(1) 

the divergenceless part of the spatial component of the four- velocity, Vj, contributes 
to the vorticity perturbation. Substituting ()3-10p and ()5-lip into ()4-100p . the com- 
ponents of the gauge invariant part ^ of the flrst-order perturbation of the 
energy-momentum tensor are obtained as 

(\^ = -i\ (5-12) 
' = (e + p) I '^u^ - (^D'^ + V* ^ I , (5-13) 

i^yr^^ = {e + p)(^D,^v' +Viy (5-14) 

(^^iJviJ. (5-15) 

Through Eqs. (f5Tl) - ([5T0|) and (f5T2]) - (f5T5]) . the hnearized Einstein equations ([2^ 
are found to be 

A^Ga^ e = i-Sndrj + A + 3K) -37i^ <P , (5-16) 
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47rGa2 (e + p) ^ A + ^ = -9^ A 1^^ -HA S +^ + 2i^) % , (5-17) 

+ 1 (^92 + 27^5^ - - ^z\^ + (ndr, + 2dr,n + n^ + $| 7, ^' 



(1) 



4a2 

+^ + 2?^a, + 2K-A) Xi' , (5-18) 



1 « 



where the component [Ti] = SirG^^^T^ * is identical to Eq. (j5-17p by virtue of 

the background Einstein equation (13-150 . 

We decompose Eqs. ()5-16P " (|5-18p similarly to Eqs. (j4-2p - ()4-4p for the metric 
perturbation hab, whose inverse relations are given by Eqs. ()4-9p ~ (j4-15p . Then, we 
obtain the equations for the scalar mode perturbations as 

inGa^ £ = {-md^ + A + ?,K) qr -SH^ <P , (5-19) 

47rGa2 (e + p) Di ^v = -5,, A 1^^ -WA (5-20) 

, (1) / , 1 \ (1) 

AnGa^ p= ( 5^ + 27^9^ -K- -Aj ^ 

+ (ndrj + 2d^n + n^ + ^z\^ (5-21) 

-A{A + m){^-^\=f), (5-22) 
the equations for the vector-mode perturbation as 



1 (1) 

AirGa^ {e+p) Vi = -{A + 2K) V, , (5-23) 

r (1) ] 

{a^{A + 2K) ^ = 0, (5-24) 



and the equation for the tensor-mode perturbation as 



[pl + 2ndr, + 2K-A) % = 0. (5-25) 



Equation (|5-25p describes the evolution of gravitational waves. 
Equation ()5-22p yields 



(1) (1) 

i^ = # . (5-26) 
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(1) 

Prom this equation, the energy density perturbation £ , the velocity perturbation 

(1) (1) 
Di V , and the pressure perturbation p are found to satisfy 

AirGa^ (Z\ - md^ - 5 {H^ - K)) (5-27) 
ATTGa"^ (e + p) A = -5^ A S -HA (5-28) 
47rGa2 V= {d^ + 3W9^ + 2dr,n + -K) ''S . (5-29) 

In the Newtonian limit, Eq. (j5-27p reduces to the usual Poisson equation for the grav- 

(1) (1) 
itational potential <P induced by the energy-density perturbation £ . This supports 

(1) 

the interpretation of ^ as the relativistic generalization of the Newtonian gravita- 
tional potential. Equation (j5-27p is the generalized form of the Poisson equation 
obtained by taking into account the expansion of the universe. 

Next, we apply the equation of state (j4-112p for the first order perturbation. 
Then, from Eqs. ()5-19p and (j5-2ip . we obtain the well-known master equation for 
the scalar mode perturbation:^)'^) 

{d^ + 3n{i + cl)dr, - clA + 2dr:n + (i + 2,cl){re - k)] 

= A-kGo^t 5 . (5-30) 

The scalar mode perturbations are completely determined by this master equation 
. (1) . 

()5-30p . If we obtain the solution <p to Eq. (j5-30p . we can obtain another scalar 

(1) (1) 

perturbation through Eq. (l5-26p . and the energy density perturbation £, the 

(1) (1) 
velocity perturbation A v , and the pressure perturbation are obtained from 

Eqs. (IF27p . (IF28I) . and (IF29p . respectively. It is also well-known that Eq. (IF30]) is 
reduced to simpler equation through a change of variables.^) 

Here, we comment on the contribution of the anisotropic stress, which is ignored 
in the above derivation of the linearized Einstein equation. If an anisotropic stress 
exists in the linear-order energy-momentum tensor, these can be formally decom- 
posed into scalar, vector and tensor types in forms similar to Eqs. (j4-3p and (j4-4p . 
If these exist, the anisotropic stress of scalar type contributes to the scalar mode 
of the perturbation and will appear on the right hand side of Eq. (j5-22p . In this 

(1) (1) 

case, the equation (|5-26p is no longer valid. Instead, — (p is proportional to the 
anisotropic stress of scalar type. As a result, the master equation (j5-30p will have a 
source term which is proportional to the anisotropic stress, in addition to the entropy 
perturbation. The anisotropic stress of vector type contributes as the source term of 

Eq. (|5-24p . if it exists at linear order. Though the solution to Eq. (I5-24P is only 
the decaying mode in the absence of the anisotropic stress of vector type, the vector 

perturbation of the metric is generated by the anisotropic stress of vector type, 
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(1) 

if it exits, and the resulting vector perturbation, Vi , of the metric directly generates 

(1) 

the divergenceless part, Vj , of the four-velocity of the fluid, which may contribute 

to the vorticity of the fluid. Finally, if the anisotropic stress of tensor type exists, it 

(1) 

appears as the source term in Eq. (j5-25p . and it generates gravitational waves, Xij ■ 

Thus, the anisotropic stress may generate all types of perturbations. The situation 

is simialr in the case of the second-order perturbations, as shown below. 

5.2. Scalar field case 

Here, we consider the linearized Einstein equation of a homogeneous isotropic 
universe filled with a single scalar field. 

We note that the background scalar field (p is homogeneous, i.e., = (p{r]), where 
T] is the conformal time. Thus, the components of the gauge invariant part of the 
first-order energy-momentum tensor '^^^ are given by 

' = (d^^idv^- S {d^ipf + a'^V'i) (5-31) 
^ = -^D^^idr^ip, (5-32) 

(^>r^ ' = ^dr,^ f ^Vi + {d^v) ] , (5-33) 



' = {dr,^idr,y.- S {d,^? - a'^V'i) • (5-34) 

Equation (l5-34p shows that there is no anisotropic stress in the energy-momentum 
tensor of the single scalar field. Then, as in the case of a perfect fluid, we obtain 
Eq. ([5^ . From Eqs. ([STD - dFTOl) . (IF26D . and (|F3T]) ~ ([5^ . the components of the 
linearized Einstein equation (j2-49p are obtained as^^ 

{+A - md^ + 4K - d^n - 271^) = 47rG {dri^id^>p:> + a^^V^i) > (5-35) 

(1) (1) 

(p +H (p = 47rG(/?i9^v9, (5-36) 

{+d^ + Sndr, + d^n + 27^2) = 4t,g (d^ipidr,f - a'^^^i^ ■ (5-37) 

In the derivation of Eqs. (l5-35p - (j5-37j) . we have used Eq. (j3-18p . We also note 
that only two of these equations are independent. Further, the vector part of the 
component ^^^G^ ^ [H] = SvrG^^^ ^ of the Einstein equations shows that there is no 

(1) 

vector mode as an initial value constraint. The equation for the tensor mode Xij is 
identical to Eq. ([F25]l . 

Combining Eqs. (j5-35p and (j5-37p . we eliminate the potential term of the scalar 
field and thereby obtain 

{d^ + A + 4K) S= 87rGdr,(pidr,ip. (5-38) 
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,f^™ (1) (1) 

Further, using Eq. (|5-36|) to express in terms of 5^ (p and <p , we also eliminate 

d-rjipi in Eq. (j5-38p . Hence, we have 

This is the master equation for the scalar mode perturbation of the cosmological per- 
turbation in universe filled with a single scalar field. It is also known that Eq. (j5-39p 
reduces to a simple equation through a change of variables, and this equation has 
the same form as Eq. (j5-30p with Cg = 1.^^ 

§6. Second-order Einstein equations 

In this section, we derive the second-order perturbation of the Einstein equation 
in the context of cosmological perturbations. In the generic case, the second order 
perturbation of the Einstein equation is given by Eq. ()2-50p . This equation is for the 
gauge invariant second-order metric perturbation Cab^ whose components are given 
in Eq. (j4-7ip . To derive the equations for each component of Cab, we have to evaluate 

In this paper, we consider the simple situation in which the first-order vector 
and tensor modes are negligible: 

= 0, Xij = 0. (6-1) 

(1) _„ 

In the linear-order perturbation, the vector mode Vi decays as a , which becomes 
smaller than the scalar perturbation. Further, as seen in N5.21 the vector mode 
(1) 

is not generated in an inflationary universe driven by a single scalar field. For 

(1) 

these reasions, it is reasonable to omit the vector mode Vi for a wide class of 

(1) 

scenarios of the evolution of the universe. By contrast, the tensor mode Xij i 
gravitational waves, may be generated by quantum fluctuations during the inflation. 
The amplitude of these stochastic gravitational waves depends on the scenario of the 
inflation. However, from the observational result of CMB,^) scalar mode fluctuations 
should be dominant, with the scalar-tensor ratio being less than unity. We can thus 

assume that the dominant contribution to the fluctuations in the universe is that 

(1) 

of scalar type. Hence, we ignore the first-order tensor mode Xij this paper. 
Of course, it is possible to extend our formulation by taking into account of the 
vector- and tensor-mode contributions to the second-order perturbations, but we 
only consider the main contribution to the second-order fluctuations in this paper. 

Because the components of the second-order gauge invariant metric perturbation 
Cab are obtained through the replacements 

(1) (2) (1) (2) (1) (2) (1) (2) ^ ^ 

Ui ^ u i, Xij Xij (6-2) 

of the variables in the first-order gauge invariant metric perturbation TLab, the eval- 
uation of ^^^^a^ [C] is accomplished through these replacements in the first-order 
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gauge invariant Einstein tensor ^^^QJ' [H], which appear in Eqs. (j5-7p - (j5-10p . Then, 
the components of '^^^QJ' [C] are given by 



(2) 



(2) 



^{-end^ + 2A + 6K) ^ -m' | , 

1 / „ (2) (2) 1 , ^ (2) \ 

[2dr^Di ^ +2'HDi <p --{A + 2K) j , 

1 f (2) (2) 1 

-2 < 2d^D' ^ +2HD' <p +- {-A + 2K + AH^ - Adr^U) 



(2) 



(6-3) 
(6-4) 

(6-5) 



1 



, , (2) (2) 

DiDJ - <P 



( (2) (2)' 

+ U2d^-A + Andr,-2K) ^ + {2nd^ + Adr^n + 2n^ + A) <p 

1 • 

+-{d^^ + 2nd^ + 2K-A) x^' ■ (6-6) 



In the simple situation descirbed by Eq. (j6-ip . the components of the quadratic 
term ^'^'^Ga^ ['^j'^] of the hnear order perturbations, which are defined by Eq. (|2-38p . 
are given by 



2 ^ (1) /(I) (1) 

-2 < I2ndr, 0/ ( - ^ 



12 I ( 



, , (1) 



(1) (1) 

8 <Ir A ^ 



m 4 r , (1) /(i) (1) , / (1) (1) , (1) 







(1) 






[2^ - 




- 2dr^Di ^ 







/(I) 




(1) i 


/(I) 




U- 


V 






-3 ^] 



(1) 



. (1) /(l) (1) , (1) 

+2DiD^ ^ [ ^ - ^] +A^ DiD^ ^ 



(6-7) 

(6-8) 
(6-9) 
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(1) . (1) (1) . (1) /(I) (1) 



(1) . (1) 



~4 A } -f^ 



(6-10) 



Through the components (j3-7p of the background Einstein tensor Gj', the compo- 
nents (j5-ip - (j5-6p of the tensor field H^^ ^, and the components (l5-7p - (|5-10p of the 
gauge invariant part ^^^Ga^ [H] of the hnear-order perturbation of the Einstein ten- 
sor, it is straightforward to confirm that ^^''Gj' [H, TC] of the components (|6-7p - (|6-10p 
satisfies the second-order perturbation of the Bianchi identity ()2-44p . i.e., 



vJ2)(?,«[7^,7^] = -2/^,,«[?^] [n] + 2h,,- [n] 

-^J^had [rL\ n Lr^ + 2Hcad [ri\ ri ■ 



(6-11) 



This imphes that we have consistently derived the components (]6-7p - (j6-10p . 

Further, we impose the relation (j5-26p . which can always be derived from the 
first-order Einstein equation when the anisotropic stress of scalar type is negligible. 
Then, the gauge invariant part of the first-order metric perturbation is given by 

Hab = -2a2 ''i^ {dr])a{dr])b - 2c? -fijidx')aidx^)b, (6-12) 
and the components (l6^- (l640]l of the tensor [n, H] are reduced as follows: 

2 



= 3D, $ +3 U 



(1) (1) 
+ 8 <P A <p 



+12 (K + W 



(6-13) 

(6-14) 
(6-15) 



2Di ^ ^ +4 <p DiD^ $ 



/ (1) i. (1) (1) (1) ( ^ (1) 
- \ 3Dk ^ D'' <p +4 ^ A <!> +[dn $ 




(6-16) 



Thus, we have evaluated '-^^^a'' [^] ^'^^ ^'^^Ga'^ [H,H] in the second-order Einstein 
equation (j2-50p . Next, we evaluate '•^^'7^ * separately in two cases, that of a universe 
filled with a single perfect fluid and that of the universe filled with a single scalar 
field. 
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6.1. Perfect fluid case 

Here, we consider the second-order perturbation of the Einstein equation in the 
case of a universe filled with a single-component perfect fluid. 

Because we concentrate only on the case in which the vector and tensor modes 
are negligible, as in Eq. (j6-ip . we should ignore the divergenceless part of the spatial 
velocity of the fluid, setting 

(1) 

V^ = 0, (6-17) 

in accordance with the first-order Einstein equation (l5-23p . Then, the components 
of the gauge invariant first-order perturbation of the fluid four-velocity are given by 

(1) (1) (1) 

Ua= -a <p {dri)a + aDi v {dx')a (6-18) 
In the simple situation that Eqs. (j6-ip and ()6-17p are satisfied, the normalization 

(2) 

condition (I4-93P of the second-order perturbation Ua of the fluid four-velocity is 
given by 

(2) (1) .(1) (2) 

u''Ua={^ \ -Di V D' V - (6-19) 

(2) 

and the components of Ua are decomposed as 

(2) / 2 (1) . (1) {2)\ / (2) (2)' 



Ua=a\^^j -D, V D' V - ,pj {dr])a + a\^D, v + V^ j {dx')a, (6-20) 
where 

. (2) 

D' Vi= 0. (6-21) 

Then, the components of the gauge invariant part *^2)7^ ^ of the second-order pertur- 
bation of the energy-momentum tensor are given by 

(2)r^ = _2(e + p)D' ^v^ A - f, (6-22) 

r9W n . U (2) (2) (1) /(I) (1) 

(2)7: " = (e + p) (a V + Vi -2 ^ A Vj + 2 f ^ + p j A V, (6-23) 

((2) (2) *^2) (1) Q\\ /(I) 

„^ -D' V -V' -6 <P D'v \ - 2(g + p] D' V, (6-24) 

(2>7: ^ = 2{e + p)Di A ''v + V % (6-25) 

Hence, through the background Einstein equations (j3-14p and (j3T5p . and the first- 
order perturbation of the Einstein equations ()5-27p - ()5-29p . the Einstein equations 



42 K. Nakamura 



for the second-order perturbations are obtained as 



(2) (2) (2) 

{-Sndr^ + A + 3K) ^ -m^ <P -AnGa^ g = Tq, (6-26) 

(2) (2) 1 ^ ^ (2) 

2drjDi W +2HDi ^ --{A + 2K) Ui 

+87TGa\e+p) (^D, %^ + V » j = li, (6-27) 

( (2) (2)1 

+ U2d^-A + mdr, -2K) <F + {2nd^ + Ad^n + 2H^ + A) ^' ^ 7i ^ 



1 ■ 

--{dl + 2Hd^ + 2K-A) Xi' 

(2) 

-SvrGa^ 7? 7.^ =r.^ (6-28) 



where 



To := ?>-KGa^{e + p)D' '^v Di %^ 

(1) . (1) / (1)^ ^ 



<p A <p -12{K + n^) l<p j , (6-29) 



9 /(I) (1) 
:= -levrGa^ f ^ + 77 j V 

(1) (1) (1) ^ (1) ^ (1) (1) 
+12n ^ Di ^ -4 ^ drjDi ^ -Adr, ^ Di ^ , (6-30) 

r. ^' := 167rGa2(e + ''^'^ ^^^^ -4A D^ S -8 D,D^ 

/ (1) . (1) (1) (1) / (1)^ ^ 
+2 SDfe ^ Z)'' <J +4 ^ Z\ ^) + I (9^ ^ 

+4 (25^7^ + ZC + 7^2W ^ j +87^^5^^ 7.^ (6-31) 

Equation (16-27P can be decomposed into scalar and vector parts. Taking the 
divergence of (j6-27p . we obtain 

8TTGa^{e+p)Di ^v= -2(9^ A -27iA S'' +DiA-^D^rk. (6-32) 
Then subtracting Eq. (j6-32p from Eq. (j6-27p . we obtain 

87rGa2(e + p) V i= + 2^) ^5 + (/^, - DiA-^D^r^) . (6-33) 
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Equation (j6-28p can be decomposed into the trace part, and the traceless part. 
This traceless part of Eq. (j6-28p can be decomposed into the scalar, vector, and 
tensor parts. The trace part of (j6-28p is given by 



1 \ (2) / , 1 \ (2) 



(2) 1 

-AirGa'^ V= -r^^- (6-34) 



The traceless part of Eq. (j6-28p is given by 



+ i + 2ndr, + 2K-A) %j = Fij - ^jijr,", (6-35) 



where Fij = 'yjkF^ ^ . Taking the divergence of Eq. ()6-35p . we obtain 

(2) (2)\ 1 ^ f 2 / ^ (2) 



-DiA + 2KDiJ - 5> j - ^dr, \a'{A + 2K) vi 

= D,rJ -\D,r^K (6-36) 
Further, taking the divergence of Eq. ()6-36p . we have 

2 /{2) (2)\ . . 1 , 

-{A^'iK)A[^ - cpj=DW,r^^--Ar,K (6-37) 

Thus, we have extracted the scalar part in traceless part (j6-35p of Eq. (|6-28p . which 
is given by 



Substituting (16-380 into ()6-36p . we obtain the vector part of Eq. ()6-28p . 



(6-38) 



dril a V. 



2 = 2a^ {A + 2K)-^ {DiA-'D'Dirf,' - D^r, '=} , (6-39) 



and substituting Eqs. (|6-38p and (j6-39p . we obtain the tensor part of Eq. (j6-28p . 
[dl + 2Hd^ + 2K-A) x^j 

= 2r,j - '^j,,r,' - 3 (^D,Dj - ^7,,/^) {A + 3K)-' (^A-^u'^DirJ - ^r,^) 

+4 (Z)(i {A + 2K)-' Dj^A-^D'OkFi - D^^ [A + 2K)-^ D'Tj)!,) ■ (6-40) 



Equations ()6-39p and (l6-40p imply that the second-order vector and tensor modes 
may be generated by the scalar-scalar mode coupling of the first-order perturbation 
if accidental cancellation in the source term does not occur. 
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Further, the equations of the scalar mode perturbations ()6-26p and (j6-34|) are 

(2) 

reduced to a single equation for (p as follows. Substituting (16-380 into (l6-26p . (l6-32p . 
and (j6-34p . the second-order perturbation of the energy density, the scalar part of 
the spatial velocity, and the pressure of the fluid are given by 

9 (2) , ^ , (2) 

47rGa2 £ = [-md^ + A + 'iK- m^) <P -Fq 



2 1 ^ ' 3 



-^Hd, {A + 3Kr' {a-^D^D.F^ ^ - ^-r,^^ , (6-41) 

(2) (2) (2) 

STTGa^{e + p)Di v = -2a,,A <P -2nD, ^ +DiA-^D^rk 

-39^ A + SKy^ (^'^D'D^r, ^ - ^rA , (6-42) 



(2) (2) 

AirGa^ V = {d^ + 3^5^ -K + 2dr,n + W^) ^ 

+3 (52 + 27^^^) + 3^)-i ^A-'DW^r^ ' - ir,'^) 

-^A-^DW,ry. (6-43) 
Through the second-order perturbation (j4-115p of the equation of state, Eqs. ()6-4ip 

(2) 

and (|6-43p yield a single equation for (p : 

(2) 

+ 37^(1 + c2)9^ - c'iA + 2dr,n + (1 + ScDin^ -K)} p 
. ^ 2 f dci ^dc'i (i)(i) dr /(i)^ 



-cUro + lrA+lU + l]A-^DW,rJ 



(a^ + (2 + 3c2) nd,) (A + 3ir)-i ^A-'DW,r^ ^ - ir,^) . (6-44) 

This is the second-order extension of the master equation (j5-30p of the scalar mode 
perturbations in the case of a universe filled with a perfect fluid. Actually, if the 
quadratic terms of the first order perturbations on the right-hand side of Eq. (l6-44p 
are absent, this equation coincides with Eq. (|5-30p . 

To solve the system of second-order perturbations of the Einstein tensor, we have 
to carry out the following process. First, we solve the first-order master equation 
(|5-30p for the perturbations. The solution to Eq. (j5-30p gives the energy density, the 
velocity, and the pressure perturbation of the fluid through Eqs. (j5-27p - (l5-29p . Next, 
we evaluate the source term of (j6-44p and solve this equation. Since the homogeneous 
solutions to Eq. (j6-44p coincide with the solutions to Eq. (|5-30p , which are known as 
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the growing and decaying modes of linear perturbation theory, the general solution 
to Eq. (j6-44p is given by an inhomogeneous solution to Eq. (j6-44p . together with 

(1) 

the growing and decaying modes of the linear-order scalar mode perturbation ^ 
with arbitrary coefficients. Once we obtain the solution to Eq. (j6-44p . we can also 
obtain the energy density, velocity, pressure perturbation at second order through 
Eqs. (|6^ - (l6^ . 

Further, we have equations for the vector and tensor modes (l6-33p . ()6-39p . and 
(j6-40p . Once we obtain the solution to the linearized Einstein equations, (j5-30p and 
(j5-27p - (|5-29p . we can evaluate the quadratic terms of the linear-order perturbations 
in Eqs. (j6-33p . (|6-39p . and (|6-40p . The evolution of the vector mode of the second- 

(2) 

order metric perturbation is determined by Eq. (j6-39p . and the rotational part V i 
of the spatial velocity of the fluid is determined by Eq. (j6-33p . The tensor mode, 
i.e., the gravitational wave mode, at second order is determined by Eq. (|6-4Up . Since 
the homogeneous solutions to Eq. (j6-40p coincide with the solutions to Eq. (j5-25p . 
the general solution to Eq. (j6-40p is also given by an inhomogeneous solution to 
Eq. (|6-40p . together with two independent solutions to (j5-25p of linear order with 
arbitrary coefficients. 

Of course, we need the additional information concerning the entropy perturba- 

(1) (2) 

tions S and S at each order to determine the first- and second-order perturbation.^^) 
Once we obtain this information, all modes of the second-order perturbation are de- 
termined by the above second-order perturbation equations, (|6-33p . (|6-38p ~ (|6-44p of 
the Einstein equation. This is one of the main results of this paper. 

6.2. Scalar field case 

In the simple situation in which the first-order vector and tensor modes are 
negligible, the components of the second-order perturbation ()4-126p of the energy- 
momentum tensor for a single scalar field are given by 



dV d'^V\ 
+Av'i^Vi + «V2^ + "'('^i)'^ ) ' (6-45) 

^% ^ = {Diip2d^ip + 2D,^idr,^i - 4D,^i drjip^ , (6-46) 



1 / (1) \ 

' = ^ I + 4a^(^ ^ Z?Vi + {drjvf +2d^ipiD'ipi ] , (6-47) 



(2Vr ^ = 2XD.;ip,D^ 
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(6-48) 



Through Eqs. (|3T8]l and (IF36]l . the components ^^^^^^ [£] + '^'^'^Qr,' [U^H] = 
SttG (2)t^ i and (^)g. [£] + (2)^. [Tt!, = SttG (2)7: of the second-order Einstein 
equation (12-50p give the single equation 

(2) (2) 1 (2) 

29^ A 'P — {A + 2K) W -SttG Di<f 2dr,f = (6-49) 



where 



:= -4(9^ <J A <^ +87^ ^ Di ^ -8 ^ dr^Di ^ +16TrGDiipid^ipi. (6-50) 



Taking the divergence of Eq. (j6-49p , we obtain the scalar part of Eq. ()6-49p , 

(2) (2) , , 

2drj +2n -8TTGip2dr,(p = A-^D^Tk- (6-51) 
Subtracting Eq. (l6-5ip from Eq. ()6-49p . we obtain the vector part of Eq. (l6-49p . 

2iA + 2K)-^ ^^DiA-^D'Tk- (6-52) 



(2) 



Through Eqs. and (ICTD . the component ^^^Grj'^iC] + ^'^^Gr^'^ [n,n] 

SttG '■2)7^ of the second-order Einstein equation (j2-50p gives 

(2) (2) 

{-Sndrj + A + 3K)^+ {-drjn - 2H^ + K) <p 



-AttG (^d^Lpdr,f2 + ^^¥^2^^ = A, (6-53) 



where 



(1) (1) (1) . (1) (1) (1) 

To := -2 ^ 52 ^ -3 a„ ^) - 3A <p D'' <p -10 ^ A <p 



-4 {dr,n + 4K + 27^2) (^^^ 

+47rG (^(a,^i)2 + Dk^iD'^^i + • (6-54) 

Similarly, through Eqs. (I3T8D and (fF38]l . the component ^^^G- ^ [/:]+(2)^. [H, H] 
8ttG ^"^^T- of the second-order Einstein equation (j2-50p yields 

/{2) (2) 

DiDj - ^ 



( (2) (2)^ 

+ I {-A + 252 + AUdr, -2K) + {2ndr^ + 2d^n + AH^ + A + 2K) <p\ 7^^- 
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+-{d^^ + 2mr, + 2K - A) Xij 



-8ttG (^dr,(pdr,(p2 " a'^(p2^{^)^ lij = ^ij, (6-55) 



where 



(1) (1) (1) (1) 
Fij := -ADi <p Dj <p -8 <P DiDj <p +167rGAv?i-Dj(^i 

+2 Z\ +4 (a^Ti + 2?^2) ^ ^ j 

+4^G (^(a^^i)2 - Dk^^D'^Vi - (6-56) 

As seen in the case of a perfect fluid in ^6.11 Eq. (|6-55p can be decomposed into the 
trace and the traceless parts. Further, the traceless part Eq. (I6-55P is decomposed 
into the scalar, vector, and tensor parts. The trace part of Eq. ()6-55ll is given by 

/ 1 \ (2) / 1 \ (2) 

( + 2718^^ --a-k \ ^ + (nd^ + d^n + 2n^ + -a + k\ <p 

-ATrG(^d^^d^^2-a^f2^{^?j = (6-57) 

where = 'y^^Fij. The traceless scalar part of Eq. (j6-55p is given by 

(2) (2) ( 1 ^ 

^ - cp=^-{A + m)-'}^A-^D'D,r^^ --r,^Y (6-58) 

where r--^:= j^^Fik. The vector part of Eq. (j6-55p is given by 

dr, (^a^ ul^ = 2a\A + 2K)-^ \^DiA~^D''DiF,^^ - D^F^ . (6-59) 
Finally, the tensor part of Eq. (16-55P is given by 

(2) 

{d^ + 2ndr, + 2K-A) Xi, 
= 2F,, - ^7.,rfc' - 3 [d,D, - i7.,Z^) + 3Kr' [a-^D^DiF^' - ^r,^) 

+4 {Z)(,(Z^ + 2K)-^D^)A-^D'DkFi ^ - D^i{A + 2RT^D''Fj)k} ■ (6-60) 
From Eqs. (j6-5ip . (j6-53p . ()6-57p . and ()6-58p . we obtain a single equation for the 

(2) 

second-order perturbation as Eq. (|6-44p in the case of a perfect fluid considered 
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in §6.11 Combining Eqs. ()6-53p and (l6-57p . we have 

/ 2 \ (2) / 1 \ (2) 

ld^^-nd.r, + -A + 2K\ ^ + indr, + -A + 2K \ 



-87rGd^ipd^ip2 = To + (6-61) 



and 



-87rGaV2|^=ro-^A^ (6-62) 



Further, substituting Eq. (j6-58p into Eq. (j6-6ip . we obtain 



(2) 

{d^ + A + AK) <p -87:Gdrjipdr,ip2 
g fc ^ * 3 

3 
2 



(6-63) 



On the other hand, differentiating Eq. (j6-5ip with respect to the conformal time rj, 
we obtain 

/ 2d?;tf „\ (2) / 2d^^ \ {2) 

- 87rGdnip2dr,^ = ( -j^drj - 2d^ j ^ + i -j^n - 2d^n - 2Hdr, \ ^ 



A-'DTk + dr^A-'D^n 



02 . 252(^ ^ \ (2) 



--^Z\-iZ)'=rfc + drjA-'D'^Tk, (6-64) 

where we have again used Eq. ()6-58p . Substituting (j6-64p into (|6-63p . we obtain the 
master equation: 

= -ro - ^A'^ + A-'DW,r^ ' + (^d, - z\-ii)'=A 

-I 1^,^ - - (Z^ + 3K)-' [A-'DW,r, ^ - ^A^} . (6-65) 
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This is the second-order extension of Eq. (l5-39p . which is the master equation of 
scalar mode of the second-order cosmological perturbation in a universe filled with 
a single scalar field. This equation would coincide with Eq. (j5-39p if the quadratic 
terms of the linear-order perturbations were absent. 

Thus, we have a set of ten equations for the second-order perturbations of a 
universe filled with a single scalar field, Eqs. (IfrKTI) . (ffr52]l . (ffr58]l - (f6^ . (ffr62]l . and 
(l6-65p . To solve this system, we have to solve the linear-order system first. Next, 
we evaluate the quadratic terms, /q) and ^ij^ of the linear-order perturbations. 
Then, using the information of the quadratic terms of the linear-order perturbation, 
we estimate the source term in Eq. ()6-65p . The general solution to Eq. ()6-65p is 
given by an inhomogeneous solution to Eq. (j6-65p and two independent homogeneous 
solutions of this equation with arbitrary constants. Since Eq. (|6-65p is the same as 
Eq. (j5-39p . except for the source term, which consists of the quadratic terms of 
the linear-order perturbations, the homogeneous solutions to (j6-65p coincide with 

(1) 

the solutions to the linear-order perturbations, ^ . Hence, we can construct the 
general solution to Eq. (j6-65p if we obtain a special solution of this equation. After 

(2) 

constructing the solution ^ to Eq. (j6-65p , we can obtain the second-order metric 

(2) 

perturbation ^ through Eq. (|6-58p . Thus, we have obtained the second-order gauge 
invariant perturbation 932 of the scalar field through Eq. (|6-5ip . Equation (j6-62p is 
then used to check the consistency of the second-order perturbation of the Klein 
Gordon equation, 

- - dV 

V"Va(^ + = 0. (6-66) 

Evaluating the source terms in Eq. (j6-60p through the evaluation of the quadratic 
terms /q, -Tj, and Fij of the linear-order perturbations, we can solve Eq. (|6-60p . We 
also note that Eq. (j6-60p is identical to Eq. (j6-40p . except for the definition of the 
quadratic terms. As in the case of a perfect fluid, this equation implies that a scalar 
mode perturbation of linear order may generat the second-order tensor mode, i.e., 
the second-order gravitational waves if accidental cancellation in the source term 
does not occur. 

(1) 

For the vector-mode perturbation ui , the situation is different from that in the 
case of a perfect fluid. Since there is no rotational component in a single scalar 
field system, there is no rotational spatial component of the velocity of the matter 
field, i.e., there is no vorticity. Instead, in addition to the evolution equation (j6-59p 
of the vector mode, we have the initial value constraint (|6-52p of the vector mode. 
However, the constraint (j6-52p and the evolution equation (j6-59p also imply that the 
second-order vector-mode perturbation may be generated by the scalar-scalar mode 
coupling of the linear order perturbations. 

Thus, we have formulated a second-order perturbation theory of a universe filled 
with a single scalar field. All modes of the second-order perturbation are deter- 
mined by the above second-order perturbation equations, ()6-5ip . ()6-52p . (j6-58p - 
(|6-60p . ()6-62p . and ()6-65p . This and the result in the case of a perfect fluid are 
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the main results of this paper. 

§7. Summary and discussions 

In summary, we have confirmed that the general formulation of the second-order 
perturbation theory developed in KN2003^^^ and KN2005^^^ is applicable to cosmo- 
logical perturbation theory. We have shown that the method for finding higher-order 
gauge invariant variables proposed in KN2003 does work in the case of cosmologi- 
cal perturbations. The key point of our method is the assumption that we already 
know the procedure for decomposing the first-order metric perturbation into gauge 
invariant and variant parts. In particular, to apply this method to higher-order per- 
turbations, we have to find the gauge variant vector field Xa in the first-order metric 
perturbation, which is defined by Eq. (j2-22p . As shown in ^4.1^ the vector field Xa 
is found by restricting the domain of the perturbations to the space of functions 
in which the Green functions (A + 2K)~^, and (Z\ + 3K)~^ can be defined, 

where A and K are the Laplacian associated with the metric jij and the curvature 
constant of the maximally symmetric three-space, respectively. As a result, we found 
gauge invariant variables for the first-order metric perturbation that are discussed in 
the literature. The resulting gauge invariant metric perturbation has the same 
form as the metric perturbation described by the longitudinal gauge (the Newtonian 
gauge). This result for linear metric perturbations was then extended to second- 
order perturbations, as proposed in KN2003. In this way, we obtain second-order 
gauge invariant metric perturbations whose components are similar to those of lin- 
ear metric perturbations in the longitudinal gauge. If we apply the gauge fixing so 
that Xa = = 0, where Xa and Ya are the first- and second-order gauge variant 
parts of the metric perturbations studied in ^ this gauge fixing corresponds to the 
Poisson gauge in the literature, which is the higher-order extension of the longitu- 
dinal gauge. Further, as proposed in KN2003, we can also define the gauge invariant 
variables for the perturbations of an arbitrary field other than the metric, as shown 
in Eqs. (l2^ and (|2^ . 

We have also shown that the formulae (|2-28p and (|2-29|) of the decomposition 
of the gauge invariant and variant parts for an arbitrary field other than the metric 
play crucial roles in the second-order gauge invariant perturbation theory. We have 
seen that the first- and second-order perturbative energy momentum-tensor of a per- 
fect fluid [Eqs. (I¥99D - (|¥T02D ] and a single scalar field [Eqs. gT23])-g426D] can be 
decomposed into gauge invariant and variant parts in the same forms as Eqs. (|2-28p 
and (j2-29p . Here again, we note that no background values of the fluid components, 
the scalar fleld, nor the metric were explicitly used in the derivation of these expres- 
sions. We thus conclude that Eqs. (I¥W|1 - (I4-102I) and Eqs. (|4-123p - (|4-126p are vahd 
in the perturbation theory on any background spacetime. 

Further, we have derived the perturbed Einstein equations of flrst and second 
order in terms of the above gauge invariant variables. As shown in KN2005, each 
order of the Einstein equations can be written in terms of only the gauge invariant 
variables, and we do not have to consider the gauge degree of freedom when we treat 
the perturbed Einstein equations. Though it is well known that the first-order met- 
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ric perturbation in a homogeneous and isotropic universe can be decomposed into 
scalar, vector, and tensor types, we have shown that the second-order metric pertur- 
bation can also be decomposed into these three types. In the perturbation theory at 
linear order, these three types of perturbations are decoupled. However, at higher 
orders, these three types of perturbations are coupled. At second order, this mode- 
mode coupling appears as a source term which consists of quadratic terms of the 
linear-order perturbations. Because the scalar mode of the perturbations yields the 
dominant contribution in many cosmological contexts, we consider only the situa- 
tion in which the first-order vector and the tensor modes are negligible. Even in this 
simple situation, the second-order Einstein equations imply that the second-order 
vector and tensor modes may be generated by the scalar-scalar mode coupling. Fur- 
ther, we have shown that the Einstein equations for the second-order perturbations 
have forms similar to those for the linear-order perturbations, but there are source 
terms due to the quadratic terms of the linear-order perturbations. Since we have 
developed the perturbation theory order by order, the Einstein equations for the 
second-order perturbations can be solved using techniques for linear-order perturba- 
tions. In particular, we have also seen that the second-order Einstein equations for 
the scalar mode perturbations are reduced to single equations in the cases of both a 
perfect fluid and a scalar field. The resulting equations have forms similar to those 
for the linear-order perturbations, but there is a source term which consists of the 
quadratic terms of the linear-order perturbations. 

Now, we discuss the definitions of the gauge invariant variables found in the 
literature. It is well known that there are many definitions of the gauge invariant 
variables for density perturbation.^^' Thus, there is no uniqueness in the defini- 
tions of the gauge invariant variables. This results from the fact that we can always 
construct new gauge invariant quantities from combinations of other gauge invariant 
variables. In many works, the interpretation of gauge invariant quantities is based 
on the coincidence of the perturbative variables in an appropriate gauge choice. For 

(1) 

example, the gauge invariant variable E defined by Eq. ^-84^ describes the den- 
sity perturbation, because the variable £ coincides with V, defined in Eq. j^- T^p in 
the gauge choice = 0. This criterion for the interpretation of a gauge invariant 
variable for density perturbations produces many different definitions of the density 
perturbations, as pointed out in the literature. ^^'^^ However, we have to recall that 
the gauge choice is the point identification map between the physical spacetime A4\ 
and the background spacetime Aio, as reviewed f|2j The concept of the gauge choice 
does not exist for the physical spacetime; it has meaning only if we introduce a refer- 
ence manifold, A4o- Moreover, because all variables on the physical spacetime, 
are pulled back to the background spacetime, Mq, these pulled-back variables nec- 
essarily depend on the gauge choice, in general. This gauge dependence is due to the 
point identification map between Aix and A^O; ^-nd it is not due to the nature A4\ 

itself. Thus, the density perturbation V in Eq. (j4-74p is defined by the pull back of 
the gauge choice and depends on the gauge choice. However, this gauge dependence 
is not due to the nature of the physical spacetime nor the background spacetime 
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themselves. We should emphasize that the physical spacetime is described only by 

the gauge invariant quantities, and hence physical density perturbation is the gauge 

. , , (1) (1) (1) 
invariant variables £ , not e . Thus, it is meaningless to interpret the variables £ in 

(1) 

terms of its coincidence with e in some gauge choice, because the physical meaning 

(1) 

of the variable e may depend on the gauge choices, whereas the physical meaning of 

(1) 

the gauge invariant variable £ does not depend on this choice. To understand why 

(1) 

the gauge invariant variable £ has the meaning of the energy density, it is enough 
to point out two facts. First, note that the gauge variant part Xa of the pulled-back 
metric perturbation hat is not a variable on the physical spacetime, which arises from 
the gauge choice, i.e., teh point identification map between the physical spacetime 
A4x and the background spacetime Mq. The information regarding the gauge choice 
is provided by the gauge variant part, Xa, of the metric perturbation, but there is no 
information regarding the physical spacetime in the gauge variant part of the metric 

(1) 

perturbation. Second, £ consists of only the variables of the energy density, its 
perturbation, and the vector field Xa of the gauge choice. Hence, the only possible 

(1) 

physical meaning of £ is that it is the energy density perturbation. Similarly, and 
more generally, the gauge invariant variables defined by Eqs. (j2-26p and (j2-27p for 
the first- and second-order perturbations of an arbitrary matter field Q have the 
physical meaning of the first- and second-order perturbation of the physical variable 
Q, respectively. 

Of course, we can ask which variable is useful when we have many definitions 
of the gauge invariant variables. This question is a different from the main point 
discussed above. To answer this question, we have to specify the problem which we 
want to clarify, and we have to specify for what the variables are useful. A partial 
answer to this question should be provided by the correspondence found between 
the gauge invariant variables and observables in observations and experiments. The 
gauge invariant variables which are useful to study physical processes should appear 
in these physical processes. Therefore, the question of which variable is useful is 
reduced to the question of which variable appears in the physical processes. Because 
observations and experiments are constructed from some physical processes, it would 
be interesting to investigate the correspondence between gauge invariant variables 
and observables in observations and experiments. Through such an investigation we 
could confirm that our formulation developed in KN2003, KN2005 and in this paper 
is relevant to physical processes. 

Ten years ago, Mukhanov et al. proposed a gauge invariant treatment of second- 
order cosmological perturbations.^^-* Their aim in investigating such perturbations 
was to clarify the back reaction effect on the expansion of the universe due to the 
inhomogeneities of the gravitational field. They also proposed an averaging proce- 
dure. In their papers, they also discussed the gauge issue in second-order general 
relativistic perturbations. In their proposal, the gauge transformation of a second- 
order perturbation should be given by an exponential map. From our understanding 
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of the "gauge" in general relativistic perturbations, which was reviewed in ^ this 
proposal corresponds to our gauge transformation with ,^(2) = 0. Moreover, in their 
treatment, the perturbative expansions of the metric and the matter fields take the 
form 

Q = Qo + A(5Q, (7-1) 

instead of Eq. ()2-10p in this paper, which includes a term of O(A^). They also 
discussed the back reaction effect due to the nonlinear effects of the Einstein equation 
through the substitution of the expansion ()7-ip into the Einstein equation, and the 
evaluation of the quadratic terms of 5Q. Thus, their treatment of the second order- 
perturbations is quite different from the formulation developed in this paper. Thus, 
we have to regard their treatment of second-order perturbations, including their 
arguments concerning the gauge, as being based on a perturbation scheme that is 
quite different from that given in this paper. Further, the correspondence between 
their works and ours are highly non-trivial. At this time, it is not clear that we should 
clarify the correspondence between their works and ours, because our formulation 
has not yet been constructed to treat the back reaction effect. If the formulations of 
the back reaction effect or some averaging procedures are formulated on the basis of 
our formulation, it will be worthwhile to compare their works and ours. 

Recently, Noh and Hwang studied second-order cosmological perturbations^^-* on 
the basis of the ADM formulation. They investigated various gauge fixing methods, 
gauge invariance, and the second-order Einstein equations in a complicated man- 
ner. Contrastingly, in our formulation, all gauge invariant variables for all fields 
were prepared before the derivation of the perturbed Einstein equations. As shown 
in KN2005, the Einstein equation is necessarily given in terms of gauge invariant 
variables only. This is shown without assuming an explicit background spacetime 
metric. Therefore, we do not have to consider the gauge degree of freedom when we 
study perturbations of the Einstein equation in both the cosmological perturbation 
theory and any other general relativistic perturbation theory, as shown in KN2005. 
In this sense, we can conclude that the formulation developed in this paper is clearer 
than the formulation of Noh and Hwang. However, it would be interesting to 
compare their approach and ours. 

In addition to the above works treating the formulation of second-order pertur- 
bation theory, there is a series of papers by Matarrese and his co-workers^^'"'^^^'^'^^ 
concerning non-Gaussian behavior generated by second-order general relativistic per- 
turbations. They also considered gauge invariant variables, but they concentrated 
on only the conserved quantities which correspond to Bardeen's parameter in the 
linear-order perturbation theory. By contast, in this paper, we found gauge invariant 
variables for the first- and second-order perturbations of all quantities. The second- 
order perturbative Einstein equations on a homogeneous isotropic background uni- 
verse were derived in terms of gauge invariant variables, without any gauge fixing. 
This is the main result of this paper. Hence, with regard to the gauge issue of 
second-order perturbation theory, we regard the formulation of second-order gauge 
invariant cosmological perturbations to be completed in this paper. Many parts of 
their works are based on the Poisson gauge explained above, and the Poisson gauge 
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is a complete gauge fixing method. Therefore, it would be interesting to follow their 
physical arguments based on the gauge invariant perturbation theory formulated in 
this paper. 

Though we have ignore the first-order vector and tensor modes, it is, in princi- 
ple, possible to include these modes. To do this, the long algebraic calculations are 
needed. Other than these long calculations, however, there is no technical problem 
to include them. If we include the first-order vector and tensor modes in our consid- 
eration, all types of mode coupling, i.e., scalar-vector, scalar-tensor, vector- vector, 
vector-tensor, and tensor-tensor mode couplings may occur, in addition to the scalar- 
scalar mode coupling discussed in this paper. These additional mode couplings are 
included in the source terms Fq, Fi, and Fij in the second-order Einstein equations, 
which consist of the quadratic terms of the first-order perturbations, as in the case 
of the scalar-scalar mode coupling studied in ^ Here, we emphasize that even if 
we consider the simple situation in which the first-order vector and tensor modes 
are negligible, the second-order vector and tensor modes may be generated by the 
scalar-scalar mode coupling if accidental cancellations in the quadratic terms of the 
linear perturbations do not occur. 

Recently, Tomita has extended his pioneering works to a universe filled with 
dust and a cosmological constant, and he also discussed non-Gaussian behavior in 
CMB due to the nonlinear effect of the gravitational perturbations. His works 
are based on the synchronous gauge. He also claimed that there is no vorticity 
perturbation in this system even if we take into account the effects of second-order 
perturbations. Contrastingly, in our analyses, the divergenceless part of the spatial 
velocity of the fluid may be generated by the non-linear effects. The divergenceless 
part of the spatial velocity of the fluid contributes to the vorticity. Of course, in 
the model of a universe filled with a single scalar field, there is no vorticity, because 
the matter current of the scalar field is proportional to the gradient of the scalar 
field. In this case, the Einstein equation gives the constraint equation for the vector- 
mode metric perturbations. Because the Einstein equations constitute a first-class 
constrained system, the initial value constraint should be consistent with the evo- 
lution equation of the vector-mode metric perturbations. In this sense, the initial 
value constraint (j6-52p for the vector mode should be consistent with the evolution 
equation (j6-59p . We can easily understand the reason for the absence of vorticity in 
a universe filled with a single scalar field, but the absence of vorticity in a universe 
filled with a dust field is not so trivial. Since the vorticity in the early universe is 
related to the generation of the magnetic field in the early universe through the Har- 
rison mechanism, and to the generation of the B-mode polarization in CMB 
anisotropy,^^) the existence of vorticity in the early universe is very important in cos- 
mology. Because vorticity perturbations of the fluid velocity related to vector-mode 
perturbations, we conclude that the generation of the vector-mode perturbations is 
an important issue in the cosmological context. 

In addition to the generation of vector-mode perturbations, the generation of 
tensor modes, which corresponds to gravitational waves, is also interesting for cos- 
mology. The upper limit of the amplitude of the vorticity perturbations and the 
gravitational wave perturbations is constrained by the observational data of CMB 
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anisotropy.'^^ However, it is known that the fluctuations of the scalar mode per- 
turbations do exist in the early universe from the anisotropy of the CMB. Hence, 
the generation of vector and tensor modes due to the scalar-scalar mode coupling in 
second-order perturbations will give a lower limits on the vorticity and gravitational 
waves in the early universe from a theoretical point of view. 

From the above discussion, we see that there are many issue which should be 
clarified using second-order cosmological perturbations. These are quite interesting 
not only from the theoretical point of view but also from the observational point of 
view. We have to clarify these issues one by one. To carry this out, the gauge in- 
variant formulation developed in this paper should provide very powerful theoretical 
tools, and we hope that these issues will be clarified in terms of the gauge invariant 
variables defined in this paper. We leave these issues as future works. 
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Appendix A 

The components of Habc [H], H^''^ [H], HJ"^ [H], and H"^" [H] 

In the derivation of the gauge invariant part of the perturbed Einstein tensor, 
the components Habc[n], H^^^[n], H^''^[n], and H"^'' [H] defined by Eqs. (I2^ 
and ()2-42p are useful. These components are given below. 
• Components of Habc [H] : 



Hr,r,rj[n] = -a%^, (A-1) 

H^nr, [n] = -a" {Di<!> + Tiui) , (A-2) 

[^] = { A<^ + {d^ + H) , (A-4) 

^ini \^ = I -^^"^^n - ^[i'^j] + ^^riXij I , (A-5) 
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• Components of [H] : 

H^vvin] = -^dr,^, (A-7) 

H,'"'[H] = -^{D,cp + nui), (A-8) 

[^] = ^ {^'^ + Wz/^} , (A-9) 

-^(9, + 2H)x/j, (A-10) 

H^^'[n] = -^{D'^+{d^ + n)u'}, (A-11) 

m = ^ [d.^lj' + I {D'l^j - DjV^) - ^9.X/} , (A-12) 

+iz?,x'' + ^t'x/'}- (A- 14) 

• Components of H^^^ [H] : 

[n] = d^0, (A-15) 

H.'^jH] = D,'P + nui, (A-16) 

iJ^^^ [H] = -D'$ - Hv\ (A-17) 

-^(9, + 2W)x,^ (A-18) 

^/^[^] = -A^-(9.; + H)zvi, (A-19) 

if^. \ [H] = dr^^^ji + D^v^^ - ^dr,Xji, (A-20) 

H^'j [H] = -dr,^^/ + ^ {D'uj - Djv^) + i5,x/, (A-21) 

Hj\ m = -^ijD^^ + 2j%Di^^ - Hj'ji^i + ^D'^xji + Dux\y (A-22) 

• Components of H"^'' [H] : 

ffvm [7^] = }^Q^^^ (A-23) 

^inv [7^] = _ 1. + ^ ( A-24) 
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(A-25) 



(A-26) 



(A-27) 



(A-28) 
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